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ABSTRACT 

We present a statistical study of several fundamental properties of radio sources in nearby clusters 
(z < 0.2), including the radial distribution within clusters, the radio luminosity function (RLF), and the 
fraction of galaxies that is radio-active (radio active fraction, RAF). The analysis is carried out for a 
sample of 573 clusters detected in the X-ray and also observed at 1.4 GHz in the NRAO VLA Sky Survey. 
The X-ray data are used to locate the cluster center and estimate cluster mass. In addition, near-IR 
data from the Two Micron All-Sky Survey are used to identify the brightest cluster galaxies (BCGs), 
and to construct the if-band luminosity function. Our main results include: (1) The surface density 
profile of radio-loud active g alactic nuclei (AGNs) is m uch more concentrated than that of all galaxies, 
and can be described by the lNavarro. Frenk. fc White] profile with concentration ~ 25. More powerful 
radio sources are more concentrated than the weaker ones. (2) A comparison of the RLFs in the clusters 
and in the field shows that the cluster AGN number density is about 5,700 times higher, corresponding 
to a factor of 6.8 higher probability of a galaxy being radio active in the cluster than in the field. We 
suggest that about 40 — 50% of radio-loud AGNs in clusters may reside in low mass galaxies (Mr > —23); 
if an equivalent population exists in the field, the RLFs in the two environments can be brought into 
better agreement. The strongest support for the existence of this low mass population comes from their 
spatial distribution and RLF. (3) The RAFs of cluster galaxies of different stellar mass are estimated 
directly. About 5% of galaxies more luminous than the characteristic luminosity (Mr- < Mk* ~ —24) 
host radio-loud AGNs. The RAF for BCGs is > 30%, and depends on the cluster mass. Compared to 
the field population, cluster galaxies have 5 — 10 times higher RAF. 

Combining the AGN RLF and spatial distribution within clusters, we estimate that they may inject an 
energy of ~ 0.13 keV per particle to the intracluster medium near the cluster center. We also investigate 
the degree of contamination by cluster radio sources on the yields of Sunyaev-Zel'dovich effect (SZE) 
cluster surveys. Under our simple model extrapolating to higher frequency, we estimate that as many 
as 10% of clusters detected at 150 GHz may host AGNs whose flux is comparable to the cluster SZE 
signal. The fraction is expected to be much higher at lower frequency. 

Subject headings: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - galaxies: active - 
galaxies: luminosity function, mass function - radio continuum: galaxies 



1. INTRODUCTION 

Studying the radio-loud active galactic nuclei (AGNs) 
is one of the primary methods to understand the central 
supermassive black hole (B H) population that is b elieved 
to be present in all galaxies (jMagorrian et al.|[l998l ). With 
the discovery that the intracluster medium (ICM) in cen- 
tral regions of clusters does not coo l to the extent pre- 
dicted by the cooling flow theory (e.g. Tamura et al.ll200lt 
iFabian et al.l 120011 : iPeterson et alj |2003[ ). more and more 
attention has been paid to the energy sources that may 
heat up the ICM. Among the interesting c andidates are the 
AGNs (in the form of radio outflows; e.g. iBinnev fc Taborl 
119951: IChurazov etal] |2001| ). Analyzing the radio prop- 
erties of the AGNs in galaxy clusters thus should pro- 
vide guidance in assess ing the fea si bility of the AGN 
feedback scenario (e.g . lEilekl 120041: iBirzan et alj [2001: 
McNamara et al.l [20051 ISanderson et alj 120051: iBest et alj 



mm wavelengths. The Sun yaev-Zerdovich effect (SZE: 
ISunvaev fc Zel'dovichl fl970h ' is a striking cluster signa- 
ture that allows clusters to be selected using high an- 
gular resolution cosmic microwave background (CMB) 
observations. In a nutshell, the hot electrons in the 
ICM inverse- Compton scatter off the CMB photons, caus- 
ing a distortion in the CMB spectrum in the direc- 
tion of a cluster. The effect is independent of clus- 
ter redshifts and therefore is an efficient way of detect- 
ing high redshift clusters (e.g. ICarlstrorn etall [20021) . 
We note, however, radio point sources that are often 
found at the centers of clusters are a potential prob- 
lem for SZE surveys. The point sources can easily be 
so powerful that they significantly contaminate the SZE 
signature from the cluster. With many SZE surveys 



2006). 



Another motivation for a comprehensive study of ra- 
dio sources in clusters is from cluster surveys at cm and 



plann ed or underway (e.g. AMI. j AMI Collaboration et all 
12001 SZA. iMuchovei et alJ l200a ACT. iKosowskvl I2003t 
SPT, iRuhl et all 12004k AMiBA. lLo et alJ l2005fT there 
are several publish ed studies that address the point 



source problem fe.g.lLoeb fc R cfregie rlll997l: ICoorav et al 



119981: iLin et~ai] [2002: Holder 20021 iMassardi fc De Zott 
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effects the cluster radio sources have on the SZE based on 
their luminosity function. 

We present here a study of several fundamental prop- 
erties of cluster radio galaxies, based on the public 
data archive of th e NRAO VLA Sky Survey (NVSS, 
iCondon et al.lll998l ). It has long been recognized that the 
two main galaxy populations that are active at radio wave- 
lengths, the galaxies with an AGN and the star forming 
galaxies, can be roughly separated by their rad io lumi- 
nosit y, with a division at ~ 10 23 WHz -1 (e.g. I Condon! 
1992). In this paper we will mainly study the properties 
of higher luminosity radio sources, the radio-loud AGNs. 
Hereafter we loosely use the term "radio source" to refer 
to the whole population, and "AGN" to refer to the high 
luminosity sources (P > 10 23 WHz -1 ). 

We study the luminosity function (LF), the radial dis- 
tribution (in terms of surface density profile), and the 
fraction of cluster galaxies that are luminous in the ra- 
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treatment, we study the radio sources within a large sam- 
ple of clusters with X-ray detections and measured red- 
shifts. In addition to ensuring the reality of the clusters, 
the X-ray data also make some important auxiliary clus- 
ter properties available, such as the cluster center and es- 
timated mass; these are used to define the cluster region 
and a fundamental length scale for the calculation of the 
radial distribution of sources. With the knowledge of the 
radial distribution, we can correct for the projection effect 
(i.e. from a cylindrical to a spherical volume) and calcu- 
late the LF. Finally, by cross-matching the NVSS source 
catalog with the source catalog of the Two Micron All-Sky 
Survey (2MASS, Uarrett et atll2000D . we probe the radio 
and near-IR (NIR) i^-band bivariate luminosity function, 
which enables estimation of the "radio active fraction" of 
cluster galaxies. These novelties and advantages of our 
study over previous ones will help theoretical modeling of 
accretion onto black holes in cluster galaxies and estimat- 
ing the selection function for SZE surveys. 

Even though both NVSS and 2MASS surveys do not 
deliver any redshift information, it is possible to construct 
ensemble cluster properties, such as the LF. Effectively, we 
are examining the properties of the excess of sources to- 
ward clusters relative to a typical location on the sky. The 
key is to know the relation between the surface number 
density and the flux (or magnitude), known as the log N- 
log S relation, derived over a large enough area so that it 
is an accurate description of the mean sky. Once we have 
this tool, we can look for cluster signals by subtracting the 
background contributions. Moreover, we can use the an- 
gular clustering properties of these sources to estimate the 
variance in this background estimate. The background- 
subtracted or residual population averaged over many clus- 
ters provides the input for our study of the cluster radio 
galaxies. This residual population is at the cluster red- 
shift. In the current study the signal is not strong because 
radio sources in clusters are not intrinsically abundant. 



One novelty of our "stacking" method over previous ones 
is that we have included the information about the volume 
surveyed for each cluster, and thus we are able to present 
the luminosity function of sources, not just the "luminos- 
ity distribution" (the mean number of sources per lumi- 
nosity interval per cluster). We have applied this method 
previous ly to study the galaxy properties in ensembles of 
clusters (iLin et al.ll20M 120041 : iLin fc Mohril2004h . 

Given the correlation between the masses of the 
central supermassive BH and the bulge of galaxies 
(Magorri an et alj Il998f) . one naturally expects that low 
mass galaxies will exhibit weaker nuclear activity. We 
present evidence suggesting the existence of a population 
of low mass cluster galaxies that host moderately pow- 
erful AGNs. By comparison of the radio LFs from all 
NVSS sources and from those that exist in both NVSS 
and 2MASS catalogs, we can infer the number density of 
AGNs that are hosted by galaxies which are too faint (not 
massive enough) to be detected by the 2MASS. 

In fj2]we give an overview of our analysis, describing the 
cluster sample selection and mass estimation, the radio 
and NIR data used, and some details necessary to carry 
out the calculations presented in later sections. The radial 
distribution of the radio sources, as well as a comparison 
with that of cluster galaxies as a whole, are presented in 
Sj3] In 2] we examine the radio and if-band LFs in de- 
tail, calculate the radio active fractions, and estimate the 
duty cycle of AGNs. In fJ5] we present evidence support- 
ing the existence of a population of low mass galaxies that 
host moderately powerful AGNs. The implications of our 
findings on the heating of the ICM and on SZE cluster 
surveys are discussed in §Sj6]and[7l respectively. Our main 
results are summarized in fJS] Unless noted, we assume 
the density parameters for the matter and the cosmologi- 
cal constant to be £Im = 0.3, — 0.7, respectively, and 
the Hubble constant to be Ho = 70/i7o km s _1 Mpc -1 . 

2. ANALYSIS OVERVIEW 

Below we first describe our cluster sample selection and 
cluster mass estimation f £|2.1[) . In £12.21 we give brief ac- 
counts of the radio and iVband data from the NVSS 
and 2MASS surveys, and the way cross-matc hing is done . 
As has been notic e d by earlier studies ( e.g. lEdgej H99ll 
ILin fc Mohn 12004 iGonzalez et~aH l2005h . the brightest 
cluster galaxies (BCGs) play an important role in clus- 
ter formation and evolution, both in the optical/near-IR 
(OIR) and radio wavelengths. We identify the BCGs in 
clusters in a statistical sense, which is described in ^2.31 
Finally in i )2.4l we examine the distribution of spectral 
index of radio sources for a subset of the sample where 
multi-frequency detections are available. This motivates 
our choice of the spectral index, which is required when 
one converts from the observed specific radio flux to spe- 
cific luminosity in the rest frame 1.4 GHz. 

2.1. Cluster Sample and Mass Estimator 

We construct our cluster sample from two large X-ray 
cluster catalogs dr awn from the ROSAT All-Sky Survey 
RASS). .NORAS (iBohringer et all 12001 and REFLEX 
Bohringer et al.ll2004 ). For brevity, hereafter we refer to 
the two samples together as the RASS sample. These cata- 
logs provide the cluster center (determined from the X-ray 
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emission) and the cluster restframe luminosity Lx ■ We se- 
lect only the clusters with galactic latitude |6| > 20 deg. 
Using the observed relatio n between the X-ray lumino sity 
and cluster binding mass (|Reiprich fc BohringeHl20020 



log 



1.4 2 L X (0.1 - 2.4 keV) 



/wf 10 40 ergs s- 



A + a log 



1AM. 



200 



(1) 

where A = —18.857 and a — 1.571, we estimate the cluster 
virial mass M200 = (47r/3)r2 o x 200p c , which is defined 
as the mass enclosed by ^200, a radius within which the 
mean overdensity is 200 times the critical density of the 
universe, p c . We note that this scal ing relation provides a 
mass estimat e accurate to 5 , 50% (IRcipric h fc Bdhringeri 
2002; see also lReiprichll2006[ ) and a virial radius r 2 oo esti- 
mate accurate to 15%. We first evolve the restframe Lx to 
z = (dividing L x by E 2 (z), where E(z) = H(z)/H and 
H(z) is the Hubble parameter), then use Eqn.[T]to obtain 
the mass and calculate the virial radius of the cluster. We 
assume any additional evolution of the scaling relation is 
negligible, as our clusters are fairly nearby (see below). 

After removing the clusters that are included in both the 
NORAS and REFLEX catalogs, those that have no red- 
shift measurements, and those that lie too close to each 
other in projection (i.e. the distance between the X-ray 
centers is less than the sum of the foci's), there are more 
than 700 clusters in the sample. We further exclude sys- 
tems at z > 0.2 or with M200 < 10 13 M Q . These cuts leave 
us with 577 clusters. 



face density er(> S). For our analysis, a relatively high 
flux threshold S u m = 10 m Jy is ch osen, which follows 
the suggestion of iBlake "fc Wall (2003)> wn0 find that near 
the NVSS completeness limit (2.5 mJy), the mean surface 
density of radio sources varies as a function of declina- 
tion. Such a fluctuation becomes insignificant at S ~ 15 
mJy. Our choice of flux threshold is a compromise between 
the uniformity of the background counts and the number 
of cluster sources available. For reference, we find that 
er(> lOm Jy) = 16.91 deg~ 2 , which is consistent with that 
found bv IBlake et all (|2004f), 16.9 deg" 2 . 
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Fig. 1. 



2.2. Radio and Near-IR Data 

To study the radio sources in clusters, we use the 1.4 
GHz data from the NVSS. The survey covers the sky north 
of —40 deg. in declination, with 45" FWHM r esolution, 
and is complete to 2.5 mJy (|Condon et al.| [l998). We use 
the survey catalog from SIMBAD to select radio sources 
within the cluster fields. Because the NVSS is carried 
out in the "snapshot" mode, sidelobes from strong sources 
may not be completely re moved during the cleaning pro- 
cess (|Condon et al.l [T998f ) . We eliminate 4 clusters with 
very strong sources (flux density S > 1000 Jy) from the 
sample. Our final sample contains 573 clusters. 

In the absence of redshift information about individual 
radio sources, we can only estimate the expected cluster 
signal statistically, which is done by subtracting the aver- 
aged background contribution from the observed data. In 
fJ3]when we study the radial distribution of radio sources, 
we determine the background level from the radial profile 
extending to large clustercentric distance. For the lumi- 
nosity function calculations (©, the background contribu- 
tion is estimated statistically from the log iV-log S 1 . With 
the statistical background subtraction, the signal-to-noise 
ratio is fairly low for most of the clusters. We therefore 
"stack" many clusters together, and we defer the descrip- 
tion of our stacking method to §[J2]and|31 

The background number is obtained by integrating the 
differential log iV-log £ above a flux threshold Su m over 
the cluster region, where the log TV-log S is derived from 
a region of 60 deg. radius centered at the north Galactic 
pole (NGP). In Fig.[T]we show the cumulative source sur- 

4 The abundance of radio stars and pulsars is very small compared to the extragalactic radio sources (Hclfand e t alJU99St) . Given that our 
clusters lie at 20 deg. away from the Galactic plane, inclusion of the point sources will increase the number of NVSS-2MASS cross-matched 
sources (by ~ 20%) without contaminating our results with Galactic objects. 



The cumulative log N-log S from all NVSS sources 
within 60 deg. from the North Galactic Pole. The NVSS survey 
completeness limit is 2.5 mJy l)Condon et a l. 1998); however, we 
analyze the sample brighter than 10 mJy throughout this paper. 

In some parts of our study we wish to connect the prop- 
erties of radio sources with the host galaxies (e.g. iden- 
tifying the brightest cluster galaxy), which in practice is 
done by cross-matching the NVSS sources with the 2MASS 
iVband point and extended source catalog^, for the K- 
band light is a good tracer of stellar mass. Preference is 
given to extended sources when there are both point and 
extended sources present within a certain angular separa- 
tion to an NVSS object. In doing so we need to deter- 
mine an optimal maximum separation d max between the 
central positions between the radio and 2MASS sources. 
We cross-correlate NVSS sources with S > 10 mJy within 
a region of 60 deg. radius centered at the NGP with all 
objects in the 2MASS catalogs and identify those whose 
positions in the two catalogs differ by less than 1'. At 
10 mJy, the p ositional uncertainty for NVSS sources is 
about 1 — 2" (|Condon et al.| [T998). Examining the dis- 
tribution of separation d for sources that are matched 
(Fig. O, we find that the distribution for NVSS-2MASS 
point source pairs shows a peak at small angular sepa- 
ration (d < 3") and a rising trend with increasing sepa- 
ration, which is due to chance projection. We therefore 
set d m ax = 5" when cross-correlating radio and iVband 
point sources. As for the NVSS-2MASS extended source 



4 



Lin & Mohr 



pairs, the distribution peaks sharply at d = 1 — 2", and 
has a roughly constant tail beyond d = 20". Because our 
clusters span a wide r ange in mass and in redshift, we 
follow the approach of Miller fc Owenl (|2001l ) and choose 
dmax for each cluster so that the probability of chance 
projections is expected to be less than 5 x 10~ 3 . In short, 
given a cluster, for every radio source within the virial 
radius, we calculate the surface density of 2MASS ex- 
tended sources around it (within 30'). The average of the 
surface densities, &k, is then used to solv e for d max via 
p(< d max) = 0.005 = 1 — exp(-7ra-_R-d^ aa .) (jCondon et all 
Il998f ). The resulting d max lies in the range o f 9 — 3 0", 
which is similar to that used bv lMiller fc Owenl (|200lh . 
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Fig. 2. — Using a maximum matching radius of 1', we first 
match NVSS sources stronger than 10 mjy with 2MASS extended 
sources. The unmatched NVSS sources are further matched with 
2MASS point sources. The resulting distributions of the separation 
are shown as the solid (NVSS with extended sources) and dashed 
(NVSS with point sources) histograms, respectively. In our analysis 
we match NVSS sources with 2MASS point sources with a maxi- 
mum separation d max = 5". See text for how d max for matching 
NVSS with 2MASS extended sources is determined. 

Our results are robust against the way cross-matching 
radius is determined; for example, using a constant d ma x — 
20" will not change the main conclusions in this paper. 

In the following, we refer to the NVSS-2MASS cross- 
correlated sources as the "XC" sources. 

2.3. Selection of Brightest Cluster Galaxies 

The most luminous cluster galaxy detected in the OIR 
waveleng ths, the so-called BCG, is a special class of object 
(see, e.g. IOemled[l97l ISchombertlll986t lDubinsklll998t 
iLin fc Mohj|2004 I Gonzalez et al. I I2005D . Their dominance 
over the cluster OIR lu minosity (especially so for low mass 
clusters, [Lin fc Mohrl 120041 ) encodes clues to the merger 
history of clusters. Outside the realm of OIR wavelengths, 
the OIR-selected BCGs are often found to be dominant 
as well (e.g. in X-ray or/and radio), and may be ac- 
tively involved in a range of thermodynamical processes 
in clusters (e.g. Tucker fc Rosnerll 1983^ Fabian et al.|[2003t 



iMcNamara et al.ll2005t ICroton et al.ll2006fh 

We identify the BCGs statistically from the 2MASS 
catalog for each cluster using two selection criteria: a 



BCG must lie within 0.1r2oo of the cluster center, and 
its iVband magnitude must be consistent with the BCG 
luminosity-cluster ma ss (L^-M) correla tion found from a 
sample of 93 BCGs (|Lin fc Mohrl [200l . We evolve the 
local Lb~M relation usin g a population synthesis model 
(|Bruzual fc Charlotll2003l) . a ssuming passive evol ution in 
luminosity for the BCG s (e.g.lStanford et al.lll998l ). As we 
have shown previously (|Lin fc Mohr]l2004[ ) , approximately 
80% of BCGs lie within 10% of the cluster center as defined 
by the peak in X-ray emission. 

With these criteria, we are able to identify the BCG 
for 342 clusters; out of these, 122 BCGs also have a ra- 
dio counterpart with S > 10 mjy. This corresponds to a 
detection rate of 36%, which is consistent with the mea- 
surements of the f raction of radio active BCGs found in 
earlier studies (e.g. lBurnsl[l990l: iBall et al.|[l99 3l). We dis- 
cuss more about the radio activity of BCGs in ij4.21 

To sec if there is any systematic difference between the 
full RASS sample and the radio-loud BCG subsample, we 
perform a 2D Kolmogorov-Smirnov test on the z v.s. Lx 
distribution for these cluster samples. The result indicates 
that the samples are similar. 

2.4. Choice of Spectral Index 

For each cluster with known redshift, we can convert 
the radio flux to luminosity. In doing so, we need to ac- 
count for the fact that what we observe at 1.4 GHz cor- 
responds to different wavelength in the restframe of the 
sources (i.e. the need to do a "^-correction" ). This re- 
quires some knowledge of the average shape of the galaxy 
spectrum in the radio wavelengths, which is typically de- 
scribed by a power-law oc v a , where v is the frequency. 

By cross-correlating the NVSS catalog with two cat- 
alogs at 4.85 GHz (from the GB6 and PM N surveys; 
IGregorv et all 119961: I Griffith fc Wright] 1 19931 and associ- 
ated papars), we can determine the distribution of spec- 
tral index a(1.4,4.85) = log(S4. 8 5/#i.4)/log(4.85/1.4) for 
radio sources within the cluster fields. For 4036 NVSS 
sources projected within the virial radius of all RASS 
clusters with S > 10 mjy, we search for their 4.85 GHz 
counterparts within 1' (note that the FWHM of the GB6 
and PMN surveys i s about 3.5' and 4.2', r espectively; 
IGregorv et all fl996l : iGriffith fc Wrightl [1991 . In total 
there are 622 identifications, corresponding to a fraction 
of 15% of all the 1.4 GHz sources considered. We caution 
that these sources may be background sources that happen 
to be projected within the cluster fields. The spectral in- 
dices are calculated using the fluxes given by the surveys; 
we do not consider effects of the differences in the beam 
size. However, we do include the uncertainties in the flux 
measurements in the determination of the mean spectral 
index. 

Because of the extraordinary properties of the BCGs, 
we also pay attention to their spectral index distribution 
(SID). Out of the 622 matches between NVSS and 4.85 
GHz cat alogs, 54 are the BCGs. Using data from the lit- 
erature (|Burnsl Il990t IBall et all 119931 ). we calculate the 
BCG SID based on 73 galaxies. 

In Fig. [3] we show the SID between 1.4 and 4.85 GHz 
for 568 non-BCG sources within cluster fields as the solid 
histogram. The SID has a mean of (5(1.4, 4.85) = — 0.47± 
0.15, with a large root-mean-square (RMS) scatter of 0.52. 
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About 17% of these sources have a > 0. The SID for the 
BCGs is shown as the dashed histogram in Fig. [3] Com- 
pared to the general population of (cluster) radio sources, 
the BCGs seem to have a steeper spectrum; the mean is 
5(1.4,4.85) = -0.74 ± 0.05, with a RMS scatter of 0.60. 
We note that 14% of the BCGs have a rising spectrum 
(a > 0) . Finally, the mean and scatter of the spectral in- 
dices from all 641 sources considered here are —0.51 ±0.15 
and 0.54, respectively. 

Because of the higher flux limits of the 4.85 GHz sur- 
veys, the SIDs deduced here may be biased against large 
negative values of a. This effect may be stronger for non- 
BCG sources as they are more abundant and potentially 
less powerful than BCGs. Partly because of this bias and 
also to re main generally consistent with pr e vious litera- 
turc (e.g. iCondonl 11992 iCoorav et all Il998l iCoble et alJ 
200(f), we adopt a = —0.8 when applying the fc-correction 
to obtain the rest frame 1.4 GHz luminosity (in §§3-6). 
In 53 however, when we extrapolate our results at 1.4 
GHz to higher frequency, we will use the SID constructed 
from all the sources considered in this section. Bearing in 
mind that a ID Kolmogorov-Smirnov test suggests that 
the probability of the SIDs of BCGs and non-BCGs being 
drawn from the same distribution is 0.003, we will discuss 
the effect on our results due to this simplifying assump- 
tion. 



cluster center. However, clusters are not of the same size 
(our clusters span a factor of ~ 200 in mass); the cluster 
mass define s a fundamental si ze scale, namely the virial 
radius fe.g. lEvrard et af] [l996). Here we examine the ra- 
dial distribution relative to r2oo, which is a proxy of the 
virial radius, and we use the X-ray luminosity to estimate 
this quantity (see Eqn.Q]). 
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Fig. 4. — Radial distribution of radio sources with luminosity 
P > 10 23 WHz _1 in clusters (solid points). The best-fit total profile 
(cluster plus background) is shown as the solid line, while the best- 
fit cluster model (an NFW profile with concentration c = 25 ± 7) is 
shown as the short dashed line. No BCGs (which, by our definition, 
all lie within 0.1r2oo) are included in the profile . Inclusion of the 
BCGs results in a profile with c = 52^^. The normalization of 
the profile is arbitrary. As a comparison we also show the surface 
density profile of all galaxies more luminous than Mx = —24 as 
the hollow points. The dotted and long dashed lines are the best-fit 
models for the total and cluster profiles. Apparently the galaxies 
have a much broader distribution than the AGNs. More details of 
the comparison are presented in £|3.3I 

The surface density of all sources within a cluster field 
can be decomposed into a cluster and a (constant) back- 
ground component. To reduce the statistical noise, we 
construct the total surface density by summing up con- 
tributions from many cluster fields. For each cluster, the 
projected distance of every source to the center (defined as 
the X-ray emission peak) is normalized to r2oo and then 
binned. We include sources out to 5r 2 oo from the cluster 
center, to better determin e the background val ue. We use 
the "universal" profile of iNavarro et al.l (|l997l . NFW) as 
the model for the spatial distribution of radio sources in 
clusters, where the normalization of the projected profile 
and the concentration c = r2oo / r s are the parameters to be 
fit. Here r s is the characteristic scale of the NFW profile. 
The best fit model and background value are obtained by 
minimizing the x 2 . The uncertainty in each radial bin is 
estimated assumi ng Poisson stat istics, following the pre- 
scription given bv lGehreld(ll986l) . 

In Fig. [4] we show the total surface density for sources 
with radio luminosity P > 10 23 WHz _1 out to 5r2oo- 
For a given lower limit of luminosity, there is a corre- 
sponding upper limit in cluster redshift z up so that all 
surface density as a function of metric distances from the the sources included have fluxes greater than our adopted 



a(1.4,4.85) 

Fig. 3. — The spectral index a for galaxies within cluster fields 
where S„ cc v a . We have matched the NVSS catalog with the GB6 
and the PMN catalogs within the RASS cluster fields. The distri- 
bution of a for 568 sources that are not associated with BCGs is 
shown as the solid histogram; the dashed histogram describes the 
SID for 73 BCGs (including 19 sources from the literature). The 
BCGs have a steeper spectral index. 



3. ENSEMBLE PROPERTIES OF CLUSTER RADIO SOURCES: 
SURFACE DENSITY PROFILE 

There have been numerous investigations of the sur- 
face density profile of the cluster rad i o sources (e.g. 
Ledlow fc Owenl Il995bt fSlee et all 119981: iMiller fc Owen 



2004; Branchesi et al 



2002t iMorrison fc Owenl 120031 : iReddv fc Yunl 12004 iBest 
2006f ). Most of these quantify the 
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NVSS completeness limit (10 mJy). The clusters used in 
constructing the surface density profile shown here are at 
z < z up = 0.0647. In total, 16,646 objects from 188 clus- 
ters contribute to the total profile, of which 836 lie within 
r2oo- At large radii, the profile is quite flat, which allows 
an accurate determination of the background level. Within 
?"200) we estimate the number of background sources is 621. 
The best-fit cluster profile (dashed line) has c — 25 ± 7. 
Mock observations are used to estimate the uncertain- 
ties in the best-fit concentration parameter and evaluate 
the systematics in the fitting procedure. We describe our 
method in QA.ll 

We note that no BCGs are included in the total profile. 
By our definition, the BCGs all lie within 0.1r2oo- Fur- 
thermore, because there are typically only a few sources 
per cluster, the BCGs stand as a relatively large propor- 
tion in the overall source sample (there are 31 BCGs that 
have log P > 23 in this 188 cluster subsample; the number 
of non-BCG cluster radio sources within r2oo is estimated 
to be 183). The inclusion of the BCGs thus makes the 
resulting profile much peakier. Fitting an NFW profile to 
it gives c = 52tu- 

3.1. Comparison with Profile from Spectroscopically 

Confirmed Members 

The surface density profile derived from our statistical 
method can be c ompared to that from a targeted spec- 
trosc opic survey (|Morrison et al.l [20031 iMorrison fc Ower] 
2003, hereafter M03a,b, respectively). M03a observe a 
sample of 30 clusters detected in the RASS out to z < 0.25 
with the VLA at 1.4 GHz in the A and C configurations. 
The sample is complete to 2 x 10 22 WHz -1 . Optical ob- 
servations are carried out for identifying and confirming 
cluster radio sources. M03b obtain the radio source sur- 
face density by stacking the clusters in the metric distance 
space, and fit the resulting profile with a King model. We 
apply our method to the data from M03a. Of the 30 clus- 
ters in the M03a sample, 14 are also in our sample. For 
these clusters we use the X-ray luminosity and X-ray cen- 
ter provided by the RASS catalogs, and construct the sur- 
face density profile based on the location of spectroscopi- 
cally confirmed members listed in Table 4 of M03a. We use 
only the galaxies with redshifts within 3000 kins" 1 from 
the cluster restframe. Within r2oo, there are 56 confirmed 
member galaxies in the ensemble cluster. The resulting 
surface density can be described by an NFW profile with 
c = 14^g 8 . Note, however, we do not distinguish and ex- 
clude the BCGs in M03's data. A fair comparison therefore 
would be to include BCGs and apply a radio luminosity 
cut at 2 x 10 22 WHz" 1 in our NVSS source sample as well. 
To reach down to P Hm = 2 x 10 22 WHz" 1 , z up = 0.0296. 
For 50 clusters that are within the redshift limit, we find 
an NFW model with c = 21^g 9 is a fair description. Fi- 
nally, if we build the surface density from the M03 sample 
with P Um = 10 23 WHz" 1 (28 objects), the best-fit has 
c = 481^2; a l so m good agreement with our result (with 
BCGs). We therefore conclude the profile derived from 
our data is consistent with that from the M03 survey. 

3.2. Dependencies on Radio Luminosity and Cluster 

Mass 



We list in Table Q] the concentration parameters of the 
best-fit NFW profiles for the radial distribution of radio 
sources selected at various luminosities, with and without 
BCGs. In the Table for fits with c > 80, we do not list 
the fit values, as these are practically power-law, and the 
NFW profile is no longer a useful description. 

From Table [1] it seems that, except for very luminous 
AGNs (those with logP > 24.5), there is no suggestion 
of luminosity-dependence of the concentration (i.e. the 
shape of the profiles). The fact that very luminous AGNs 
are more centrally located has long been noticed (e.g. 
Ledlow fc Owenlll995bt iMiller fe Ower]|2002t iRizza et all 
2003J, M03b)7~ This can result from the combined effects 
of (1) more luminous galaxies are more strongly clustered, 
and (2) more luminous galaxies have larger probability of 
hosting more powerful AGNs. We will discuss these in 
H3.31 A third possibility is the confining pressure due to 
the dense ICM near the cluster center. In H5.4I we will 
discuss this issue in more detail. 

Earlier studies of the radio source surface density have 
usually compared the data to a King profile, and found 
evidence of an excess of sources in the central bins. This 
suggests that the distribution is more concentrated than 
the model. As we have shown here, a cuspy profile like 
the NFW (or even a power-law) provides a reasonable de- 
scription of the composite surface density profile of cluster 
radio galaxies. 

We also check if there are any cluster mass-related 
trends in the radial distribution. We consider two limiting 
radio luminosities to examine the issue (logP/; m = 23.5 
and 24), and arbitrarily set the division between the high 
and low mass clusters at logMaoo = 14-2. The result- 
ing best-fit concentration parameters are listed in Tabled 
When BCGs are included, the profiles are very different 
for low and high mass clusters, although the difference is 
much smaller when BCGs are excluded from the fit. This 
probably mainly reflects the different radio active fraction 
for BCGs in high and low mass clusters (i.e. BCGs in high 
mass clusters tend to be more active in the radio; see H4.2[) . 
We note in passing that, with log Pu m — 24, we can in- 
clude clusters with redshift up to z up = 0.1899, which 
allows us to check if there is any redshift-dependent trend. 
Separating the sample at z = 0.1, we find that the radial 
distributions at the two redshift bins are similar. 

3.3. Comparison with K-band Galaxy Radial Profile 

Following th e same procedure outlined in |J3] (also see 
iLin et alj|200l hereafter LMS04), we find the best-fit clus- 
ter profile from the total surface density for 2MASS galax- 
ies. T he fc-correction of the form k(z) = — 61og(l + z) is 
used (jKochanek et al.l l200lt ). The galaxies are selected 
from a cluster subsample which is constructed so that no 
cluster lies in projection within six times the virial radius 
of any other cluster. We refer to these clusters as the "iso- 
lated" sample, and will explain below the reason for using 
this sample to study the profile. A comparison between 
the galaxy total surface density profile (Mk < — 24; hol- 
low points) with that of the radio sources (logP > 23; solid 
points) is shown in Fig. 2] The best-fit NFW profile has 
c = 4.2lo|- The normalization of the profiles is arbitrary, 
and we have shifted the radio source profile so that the 
best-fit cluster models (long and short dashed lines) meet 
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Table 1 

Radio Luminosity Dependence of Radial Distribution 



concentration 





z b 

up 




w/o BCGs 


w/ BCGs 


22.5 


0.03705 


84 


25 ± 12 


*°-15 


23.0 


0.06470 


188 


25 ±7 


r 9 +25 
oz -14 


23.5 


0.1118 


364 


21 ±7 


48±il 


24.0 


0.1899 


554 


22 ±6 


66±?j 


24.5 


0.20 


573 


59 ±11 


> 80 


25.0 


0.20 


573 


> 80 


> 80 



Note.- Surface density constructed from NVSS-only sources, using the full RASS 

cluster sample. 

a Minimum radio luminosity. 

b Maximum redshift for a cluster to be included in each subsample. 
c Number of clusters. 



Table 2 

Cluster Mass Dependence of Radial Distribution 









concentration 


log Pum 


log Af 20 o 


N c i 


w/o BCGs 


w/ BCGs 


23.5 


13-14.2 


106 


151- 


17 +1 7 ° 




14.2 - 15.5 


258 


23^ 


61+28 


24.0 


13- 14.2 


106 


13±i' 


14jf 




14.2- 15.5 


448 


24 -5 


79±lf 



Note.- Surface density constructed from NVSS-only sources, using the full RASS 
cluster sample. 



in the innermost point in the figure. It is apparent from 
the hgure that the two populations have very different ra- 
dial distributions. We note that by fitting a King model 
to the projected surface density profiles, M03b find that 
the radio-loud AGNs (HRLGs in their terminology) has a 
smaller core radius compared to the red galaxies, implying 
that powerful radio galaxies is more centrally distributed. 

Similar to our treatment of radio sources, we use mock 
observations to estimate the uncertainties in the best-fit 
concentration parameter; the details of the procedure are 
provided in ^A.2l 

Compared to the galaxy concentration c — 2.9 ± 0.2 
found in LMS04 (§3.1 therein), the value obtained here 
(c = 4.2l ; 4 ) seems to be quite high. This discrepancy is 
caused by three factors: (1) the use of a luminosity cut in 
selecting the galaxies in the present analysis (as opposed to 
including all galaxies brighter than the 2MASS complete- 
ness limit), (2) the use of a different cluster mass indicator 
(X-ray luminosity as opposed to X-ray temperature), and 
(3) the sample selection and background estimation. We 
explain their effects in turn. 

i. Luminosity Cut: As we will show below, more luminous 
galaxies are more centrally distributed in clusters. There- 
fore, mixing faint and luminous galaxies effectively reduces 
the concentration of the surface density profile. Further- 
more, as the faint galaxies are more numerous than the 
luminous ones, they have more weight in the resulting pro- 



file. 

ii. Cluster Mass Estimators: It has been found that 
for clusters hotter than ~ 5 keV, the mass derived from 
Lx — M relation is about 30% larger th an that from the 
T x - M relation (jPopesso et all |2005[ ). For the same 
galaxy distribution (in metric space), a larger virial radius 
(^200 ) corresponds to a more centrally concentrated pro- 
file, namely a larger c. With the difference in mass caused 
by the different scaling relations, we estimate the change 
in concentration is at 10 — 20% level (see e.g. §6 of LMS04 
for discussion on the effects of cluster mass uncertainty in 
the estimation of concentration). To further quantify this 
effect, we compile a sample of 50 clusters which are present 
in both the RASS sample and the 93 cluster sample used 
in LMS04 (and therefore both Lx and Tx are available). 
We construct two surface density profiles with this sample, 
one using Lx to infer 7-200 and the other using Tx- It is 
found that the concentration value of the profile using Lx 
as the mass indicator is higher than its Tx counterpart by 
15%, in good agreement with the estimate above. 
in. Cluster Sample Selection and Background Estimation: 
Another difference of the two analyses is the way back- 
ground level is estimated. In LMS04 the number of back- 
ground galaxies is estimated based on the global 2MASS 
number count (log iV-log S) . In the current analysis, we 
fit the cluster and background surface densities simultane- 
ously, using galaxies extending to 5r2oo- While the former 
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method ignores the presence of large scale structure in the 
foreground or background of the clusters, the latter (which 
we refer to as the "local" estimate) may be affected by it 
strongly. Averaging over many different line-of-sights, one 
might hope to find an agreement between the two esti- 
mates. In §3TT] we show that this is indeed the case for 
the radio sources with the RASS cluster sample (which 
is the reason we use the full cluster sample to study the 
profile of the radio sources in . Most likely because of 
their extended spatial distribution, for galaxies, the com- 
parison of the global and local background estimates shows 
a surprising difference. For Mk < —24 galaxies selected 
from the full RASS cluster sample, we find that within 
three, four, and five times the virial radius, the number 
of background sources derived locally is 14%, 23%, and 
24% higher than that based on the global estimates. To 
better understand the origin of the discrepancy, we have 
selected a subsample from all the RASS clusters which are 
not closer to each other (in projection) than 8r2oo, and 
found that the differences between the locally and glob- 
ally derived background estimates become < 4% beyond 
3r2oo- Based on the expectation that the two methods 
of background estimation should converge at some large 
clustercentric distance, we regard such a cluster sample as 
suitable for the study of the surface density profile of clus- 
ter galaxies. The drawback of demanding clusters to be 
isolated out to 8r2oo is that the sample size becomes small 
(~ 50) and the measurement of the profile is noisy. Wc find 
that a subsample isolated out to 6r2oo is a good compro- 
mise between the accuracy of background determination 
(< 5%) and the number of clusters available (~ 100, for 
studying Mk < —24 galaxies). 

In short, we attribute the differences in the value of con- 
centration found in the current analysis and our previous 
study to the three factors described above. With a better 
understanding of the systematics of the procedure to fit the 
surface density profile, the results presented here should be 
more representative of the nature of cluster galaxies. 

To conclude this section, we examine the dependen- 
cies of the galaxy concentration on the magnitudes of the 
galaxies. The results are shown in Table [3] The last row in 
Table [3] suggest a trend of luminosity segregation; the most 
luminous galaxies (e.g. Mk < —25) tend to distribute with 
larger concentrations. In LMS04 (§5.2 therein) we examine 
a sample of 93 clusters and find that there is no difference 
in concentration for luminous or fainter galaxies (with an 
arbitrary division set at Mk = —23.5). We repeat the 
analysis by examining the concentration for galaxies with 
M K < -23.5 and -23.5 < M K < -21. We find that 
c = 4.2+fj;4 and 3.9i ;|, respectively. Essentially, the lumi- 
nosity segregation is only strongly apparent with the most 
luminous galaxies, and the larger cluster sample used here 
allows analysis of more finely binned gradations in galaxy 
luminosity. 

4. ENSEMBLE PROPERTIES OF CLUSTER RADIO SOURCES: 
RADIO AND K-BAND BIVARIATE LUMINOSITY 
FUNCTION 

Consider the K-band and radio bivariate luminosity 
function of galaxies tp, where i/j(Lk, P)dLK<iP is the spa- 

5 Contaminations due to point sources also happen in X-ray surveys, 
with s pectroscopic confirmation, a s well as correlations with galaxy 
small iBohringer et alj|2000t [200l) . 



tial number density of galaxies having if-band and ra- 
dio luminosities within the intervals [Lk, Lk + ALk\ and 
[P, P + dP), respectively. When this function is projected 
onto one of the luminosity axes, one obtains a univariate 
LF 4>{L K ) or 4>{P). 

The bivariate LF %1){Lki P) is extremely useful, because 
it can be thought of as the conditional probability of a 
galaxy of a given .fT-band luminosity (a mass surrogate) 
to host an AGN of radio luminosity P. With it we can be- 
gin to explore how the central engines are fueled (whether 
galaxy environment is relevant?) and how long outbursts 
or periods of activity last (i.e. the "duty cycle"), which in 
turn sho uld provide important insights into galaxy evolu- 
tion (e.g. ILedlow fc Owen|[l996l hereafter L096). 

Without accurate redshift determinations, it is not pos- 
sible to map out the bivariate LF (BLF) fully. Using the 
NVSS-2MASS XC sources, however, helps significantly as 
they allow us to navigate through the Lk~P space. By 
specifying a cutoff in cither the radio or if-band fluxes, we 
can build semi-BLFs such as iP(Lk, > P) and ip(> Lk, P)- 
Furthermore, comparisons between these semi-BLFs and 
the univariate LFs yield clues about the populations that 
have either fT-band or radio luminosities too faint to be 
included in the surveys. For example, the ratio between 
4>{Lk,> P) and 4>(Lk) is the fraction of galaxies with 
a given Jf-band luminosity that have radio luminosity 
greater than P, which we refer to as the radio active frac- 
tion (RAF). 

The pioneering study of L096 (see also lLedlow fc Owenl 
19953B) examines the BLF based on a large sample of 
Abell clusters. Our results build upo n their results for sev- 
eral reasons. First. ILedlow fc Owenl ' s analysis, despite its 
large cluster sample and the spectroscopic redshift confir- 
mation of the cluster memberships, focuses only on sources 
within 0.3 Abell radius, and thus may not reflect the full 
cluster population. Second, they adopt the cluster optical 
LF from other studies, while we construct the JC-band LF 
from the same sets of clusters. Third, our cluster sample is 
X-ray selected, thus significantly reducing contamination 
from spurious clusters that are due to chance projection^. 
In addition, X-ray data provide an estimate of the cluster 
virial region, which sets a physically motivated region for 
the examination of cluster properties. Finally, we pay par- 
ticular attention to any possible trends with cluster mass, 
and to the effects of the BCGs, neither of which have been 
addressed in previous studies. 

4.1. Radio Luminosity Function 

We construct the radio source LF based on a method 
developed in LMS04. Essentially, we use the cluster radio 
galaxy excess toward each cluster to build up a luminosity 
function. Specifically, for each cluster, we assume every 
source is at the cluster redshift z c and transform the flux 
density into the luminosity (in units of WHz -1 ) with the 
standard fc-correction (e.g. M03a), assuming an average 
spectral index of a = —0.8 (see S}2T]). We bin the radio 
sources projected within r2oo in the logarithmic luminosity 
space. The faintest bin included is the bin with a faint edge 
that is just greater than the luminosity that the limiting 
flux density (i.e. 10 mJy) corresponds to at z c . For each 

However, as the NORAS and REFLEX surveys have been complemented 
distribution, the contamination rate in the final catalogs should be very 
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Table 3 

Tf-BAND Luminosity Dependence of Radial Distribution 



concentration 



MK,Um 




N cl 


w/o BCGs 


w/ BCGs 


-22.0 


0.02710 


15 


5 +i - 3 

°- u -1.0 




-23.0 


0.04325 


49 


o fi +0.5 
o.D_ 4 




-24.0 


0.06930 


106 


4.2t° j 


6 4 +0 - 7 


-25.0 


0.1117 


207 


7.8±{1 


28.511;? 



Note.- Surface density constructed from 2MASS sources, using the "isolated" cluster 
subsample (see E)3.3|l . 



bin we estimate the background count from the logiV- 
log S (by transforming the luminosity interval into flux 
densities). Once every cluster field has been examined, for 
each bin we subtract the cumulative background contribu- 
tion from the total observed number of sources, and that 
difference is then divided by the sum of volumes from all 
the clusters that contribute to the bin. For each cluster we 
adjust the cluster volume to account for the expected dif- 
ference between the number of cluster sources identified in 
the cylindrical volume surveyed and the number expected 
in the spherical volume we are trying to study. We do this 
using an NFW profile for the cluster radio galaxies with 
c = 25, consistent with the observed radial profile (see 
e.g. Table [1}. 

The uncertainty in each bin includes the Poisson noise, 
the cosmic variance in the background count, and the un- 
certainty in cluster mass. We foll ow the standard treat- 
ment fe.g. lWhite k. Srednickilll995l ) for the variance in the 
mean surface density of the background 

167H(1 - cos6»2oo) 2 ^ 
Bi=2ir f Pi{x)dx, (2) 

J COS 9 2 00 

where a is the mean background surface density (from 
the log iV-log S) , 0200 is the cluster angular si ze, C; is the 
radio source angular power spectrum (from iBlake et alj 
2004), Pi(x) is the Legendre polynomial, and l m ax is the 
maximum angular frequency corresponding to the angular 
resolution of the NVSS. The cosmic variance in each lumi- 
nosity bin is calculated by summing up in quadrature the 
a cv terms from every cluster. 

To account for the effect of uncertainties in cluster mass, 
given a cluster sample, for every cluster we perturb its 
mass from the nominal value (inferred from Lx) in the 
logarithmic space by a Gaussian random number with a 
standard deviation of 0.2 (appropriate for a fractional un- 
certainty in mass of ~ 50%; tJ2.1[) . and use the resulting 
?*200 to include radio sources and construct the LF. Re- 
peating this process ten times, for each luminosity bin, 
we calculate the dispersion of the ten LFs, which is then 
summed with the uncertainties due to Poisson statistics 
and the cosmic variance. 

We show in Fig. [5] the composite radio LF (RLF) for 
all 573 clusters (solid points). To conveniently describe 
our data, we take the LFs for AGN and star forming 



(SF) galaxies determined by ICondon et al.l (|2002t here- 
after C02) for the field population and fit the sum of the 
two to the data, allowing the amplitudes of the RLFs to 
vary, while fixing the shape of the RLFs. In doing so we 
are assuming that the shape of the RLFs are the same 
both within and outside clusters. The reason we only fit 
for the amplitudes is because our data only allow us to 
fit for the combined contribution from the AGN and SF 
populations. The functional form that C02 adopt is 

\og<P = y-^+ ^° gP w ~ X ^ ' -l.SlogP. (3) 

For the AGN RLF, C02 find that (b,x,w) = 
(2.4, 25.8, 0.78); the same set of parameters for the SF RLF 
is (b, x, w) = (1.9,22.35,0.67). These parameters, which 
control the shape of the RLFs, are kept fixed during fit- 
ting. We find that RLFs with the amplitude parameters 
y = 37.97 and 35.00 best describe the cluster AGN and 
star forming galaxy populations, respectively. The dotted 
line in Fig. [5] shows the sum of our best-fit cluster AGN 
and SF RLFs. 

In Fig. we also plot the field RLF from C02 (AGN 
and SF galaxies combined), multiplied by a factor of 

200/fi M (~ 855 ±ioi> where 

we use the value favored 

by the three-year WMAP result, Q M = 0.238±g;^, 
Spergel et al.1 I2006D to account for the difference in the 
mean overdensities in the clusters and the field (dashed 
line). Compared to the (scaled) field RLF, there is an over- 
abundance of radio sources in clusters. As a way to quan- 
tify the difference, we integrate both the cluster and (un- 
sealed) field AGN RLFs over the same range in the radio 
luminosity. Within the interval log P = 24 — 27, the ratio 
between the resulting number densities is « 5690 ± 1000, 
or 6.8 ± 1.7 time s the differen c e in o verdensity (200/Om)- 
Using the RLF of lSadler et ail (|2002t) results in a larger dif- 
ference (w 8690 ± 1520). In these calculations we have as- 
sumed a 5% uncertainty of th e field RLFs. Our RLF iscon- 
sistent with that obtained bv lMassardi fc De Zottil (|2004[ ); 
integrating their RLF (as tabulated in Table 2 therein), we 
find the number density in clusters is w 5540 ± 1090 times 
higher than that in the field (corresponding to 6.6 ± 1.8 
times the difference in overdensity). 

The comparison is robust against our choice of the mean 
spectral index. Using a — —0.5, which is the mean derived 
from all the sources for which we have both 1.4 and 4.85 
GHz fluxes ( N2.4p . we find the ratio of number densities is 
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6.7 ± 1.7 times higher than 200/O 



M • 




22 24 26 

log[P 1 .4GH Z /( W Hz" 1 )] 

Fig. 5. — The composite radio luminosity function (RLF) within 
7*200- We denote "per log 10 interval of radio luminosity" by "dex — 1 ". 
We construct the RLF based on two data sets: the NVSS (solid 
points) and FIRST (hollow points) surveys. The two RLFs agree 
with each other (see JfTTJl. The dotted line is the fit to the NVSS 
RLF (see text for m ore details). The dashed line is the field RLF 
llCondon et al-lfeOO^ , shifted in amplitude to account for different 
overdensities in the field and cluster environments. There is an over- 
abundance of radio sources in clusters, compared to the field envi- 
ronment. We do not show the lower-side errorbars for bins whose 
ler value is consistent with zero. 

Our finding that radio galaxies are over-abundant in 
clusters is at odds with the conclusion of L096, who find 
no differences between the fraction of all galaxies that is 
active in the radio in and outside clusters. While our re- 
sult is based on the sources solely identified in the ra- 
dio (irrespective of their optical property), in the analysis 
of L096, there is an explicit requirement that the radio 
sources must have an optical counterpart. This is one of 
the main reasons that lead to the different conclusions. 
A fair comparison to L096's analysis, therefore, may be 
done by using the sources that exist in both NVSS and 
2MASS catalogs (our "XC" sources; We will com- 

pare our results with theirs and comment on this issue in 
later sections (see H4.4p . 

Before concluding the section, we check for any cluster 
mass dependencies of the RLF. We separate the clusters 
into high (log M 200 > 14.2) and low mass samples. The 
most luminous AGNs (logP > 25.5) are only found in high 
mass clusters. On the other hand, for less luminous bins, 
the amplitude of the RLF for low mass clusters seems to 
be slightly larger. Summing up the RLF for luminosity 
bins with logP > 23 gives the number density of AGNs; 
we find that the AGN number density in low mass clusters 
is 2.0 ± 0.7 times higher than that in high mass ones, that 
is, the data suggest that low mass clusters has a higher 
amplitude of the RLF. 

4.2. K-band Luminosity Function and Radio Active 
Fraction 

6 The cluster sample used in this section is the full RASS sample; we 
essentially unchanged. 



Here we study the /f-band LF (KLF) of the radio 
sources, and compare it to that from all 2MASS galaxies. 
The statistical background subtraction method is similar 
in spirit to the one adopted in the previous section; see 
LMS04 for more details. This analysis will provide in- 
sights into the probability of a galaxy hosting a radio-loud 
AGN. 

We can only obtain KLF for NVSS sources that have 
counterparts in the 2MASS catalog, namely the XC 
sources. We show in Fig. [5] the KLFs for 2MASS and 
XC sources, constructed from 364 cluster^ at z < 0.1118. 
The upper limit of redshift is set so that all the XC sources 
have logP > 23.5. The solid symbols are the 2MASS- 
only sources, and the hollow points are the XC sources. 
To correct for the background in the XC sources, at each 
cluster redshift we calculate the flux density Si that cor- 
responds to the limiting radio luminosity, and use the K- 
band log iV-log S constructed from the XC sources whose 
radio flux density S > Si. 



u 
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Fig. 6. — The composite K-band luminosity function within r2oo 
from 364 clusters at z < 0.1118 (solid points). The cross-correlated 
sources (with logP > 23.5) are sh own as hollow symbo ls. Also 
shown as the solid line is the field LF (Kochanek ct al. 2001), shifted 
vertically to account for different overdensities in the field and clus- 
ter environments. We do not plot the errorbars for points whose 
1-cr uncertainty is consistent with zero. The inclusion of BCGs in 
the LF makes it deviate from the Schechter function at the bright 
end. 



Also shown is the field LF from lKochanek et alJ (|200ll ). 
scaled up by 200/17 m to account for different overdensi- 
ties in the field and cluster environments. Using the KLF 
determined in LMS04, which is a good description of the 
data here, we find the cluster-to-scaled field number den- 
sity ratio is 1.02 ± 0.22 down to M K = -21. 

By integrating the 2MASS-only and XC KLFs to a given 
K-band magnitude Mn m , we can calculate the fraction 
of galaxies more luminous than Mu m that are hosts of 
P > Piim radio sources. As is clear from the Figure, 
the radio active fraction (RAF) is a strong function of 
K-band luminosity or galactic stellar mass. We find that 
about 11.7% of massive galaxies {Mk < —25) have ra- 

have checked that using the "isolated" subsamplc (sec i|3,3H the KLF is 
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dio luminosity logP > 23.5, but only 2.8% (0.74%) of 
galaxies more luminous than —24 (—22) mag host AGNs 
of the same radio luminosity range (logP > 23.5). Ta- 
ble [U contains the RAFs for 22 < log Pu m < 25 and 
-22 > M K ,Um > -26. In the Table we include RAFs cal- 
culated when the BCGs are included or excluded (shown 
in parenthesis). Note that at a given Pu m , a correspond- 
ing maximum rcdshift z up defines the cluster subsamplc 
used to calculate the RAF, and therefore different entries 
within a Mkmiu column are not calculated based on the 
same subsample. 

A couple of trends are apparent in Table [D First, for 
galaxies of the same mass class (i.e. same Mk Um)i fewer 
galaxies are radio-active when Pu m increases. Second, the 
cumulative RAF increases dramatically for a given Pu m 
when considering more if-band luminous (i.e. more mas- 
sive) galaxies. For example, 4.9% of the galaxies of similar 
or higher luminosity than M* host AGNs (log P > 23). For 
galaxies with Mk < —26, the fraction becomes 45.3%. If 
we exclude the BCGs from the calculation, the luminous, 
non-BCG galaxies still have « 5 times higher probability 
(RAF = 16.4%, compared to the 3.3% RAF for M K < -24 
non-BCG galaxies) to host radio-loud AGNs. 

We can calculate the RAF specifically for the BCGs, 
as the ratio between the number of BCGs having radio 
luminosity P > Pu m and all the BCGs that can be identi- 
fied statistically (342) in the RASS cluster sample. For 
logP /im = 23, 24, 25, and 26, we calculate the BCG 
RAFs and record them in Table [5] Compared to non- 
BCG galaxies of M K < M*, the BCGs are w 10 times 
more likely to host an AGN with logP > 23, and ~ 29 
times more probable to harbor logP > 25 AGNs! To find 
out how much more likely it is for BCGs to host radio- 
loud AGNs compared to galaxies of comparable luminos- 
ity/mass, we compare the RAFs for BCGs more luminous 
than Mk — —26 (Table to those listed in parenthe- 
ses under the M K < -26 column in Table H The BCGs 
are about 2 — 3.6 times more likely to be active in the 
radio than the similarly massive, non-BCG cluster galax- 
ie s. These findings are in agreement with those presented 
in iBest et ail (|2006l ) . This boost in radio activity is most 
likely due to the special location the BCGs have compared 
to other galaxies, namely at the center of the cluster. 

We note in passing that our results are not affected 
by the choice of the spectral index (a = —0.8); using 
a = —0.5 lowers the RAF very slightly. For the most 
powerful sources (e.g. logP > 25) for which the effect is 
strongest, the change in RAF is only at few percent level. 

4.2.1. Dependencies on Cluster Mass 

Next we examine the RAF as a function of cluster mass 
(when the limiting radio and if -band luminosities are kept 
the same), both for the whole galaxy population and for 
the BCGs. In Fig. [7] we show the comparison of RAFs for 
two limiting radio luminosities (log Pu m — 23, open sym- 
bols, and log Pi m = 24, solid symbols) as a function of 
limiting if -band magnitude MK,Um- The division of high 
and low mass is arbitrarily set to be logM2oo = 14.2 The 
RAFs of high (low) mass clusters are shown as circles (tri- 
angles). For MK,Um > —24, the RAFs of high and low 
mass clusters are similar; for galaxies more luminous than 
Mk = —24, those in high mass clusters appear to be more 



active in the radio than their counterparts in low mass 
clusters. While the RAFs in high and low mass cluster 
differ at < 2a level for the case of log Pu m = 23, the differ- 
ences become more significant for the case of log Pu m = 24 
(> 3c; e.g. compare filled circles with filled triangles). 
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Fig. 7. — Comparison of radio active fractions (RAFs) in high and 
low mass clusters, where the division in mass is set at log M200 = 
14.2. The RAF is denned as the fraction of galaxies more lumi- 
nous than Mji n m which also host radio sources more powerful than 
Plim- The RAF as a function of n m of high (low) mass clusters 
is shown as circles (triangles). The open and solid symbols denote 
the cases of log Pu m = 23 and 24, respectively. For both cases the 
high mass clusters have higher RAFs when Mj^ i irn < —24. 

That RAF weakly correlates with cluster mass seems 
to be counter-intuitive to the finding that the RLF of the 
low mass clusters has a slightly higher amplitude than that 
of the high mass clusters (347TJ) . This is because (1) for 
massive clusters the abundances of both luminous galax- 
ies and radio-loud AGNs at the bright-end of the KLF are 
higher than those in low mass clusters, and (2) because of 
the shape and amplitude of A"-band and radio LFs, there 
are many more galaxies (Mk Mk*) than AGNs in low 
mass clusters. Interestingly, these reflect the main differ- 
ences between the KLFs in high and low mass clusters (the 
lack of very luminous galaxies and the higher amplitude 
of the KLF near M* in low mass clusters; see LMS04). 

The trend for a higher RAF in more massive clusters is 
stronger for the BCGs. Out of the 342 clusters for which 
we can identify a BCG statistically, 288 of them are more 
massive than logM2oo = 14.2, and 54 are lower mass clus- 
ters. Among the BCGs in high mass clusters, 36 ± 3% 
have radio luminosity logP > 23, while only 13 ± 5% of 
the BCGs in low mass clusters are as powerful in the radio. 
The average radio luminosity for the BCGs in high mass 
clusters is much higher than that for the low mass clus- 
ter BCGs (geometric means are 10 24 - 30±0 07 WHz~ 1 and 
10 22.88±o.i8 WHz -i ; respectively). This enhanced RAF for 
BCGs in more massive clusters may provide a way to assess 
the relative importance of gas supply, local ICM pressure, 
and interactions/mergers between galaxies in clusters of 
different masses. 

Finally, we compare our results with the findings of lBestl 
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Table 4 
Radio Active Fraction (%) 



log Pu m 


Mk < -22 


M K < 


-23 


M K < -24 


M K < -25 


M K < -26 


22 


8.1 (7.2) 


12.5 


(11.0) 


19.9 (16.7) 


45.7 (35.7) 


37.6 (18.3) 


23 


1.5 (1.1) 


2.1 


( 1-4) 


4.9 ( 3.3) 


19.1 (11.6) 


45.3 (16.4) 


24 


0.42 (0.26) 


0.66 


( 0.39) 


1.6 ( 0.95) 


7.5 ( 4.0) 


19.1 ( 5.9) 


25 


0.093 (0.048) 


0.16 


( 0.082) 


0.40 ( 0.20) 


2.0 ( 1.0) 


6.6 ( 2.3) 



Note.- All values calculated within 7-200- Values in parentheses are calculated without BCGs. See Table [5] for BCG-only RAFs. 



Table 5 

BCG Radio Active Fraction (%) 



log P Um 


all 


M K < -26 


23 


32.7 


35.2 


24 


19.9 


21.0 


25 


5.84 


6.85 


26 


0.585 


0.913 



(2003), who suggests that the radio-loud AGNs are prefer- 
entially located in groups and lo w ma ss clusters. By cross- 
matching NVSS with 2dFGRS, IBestJ examines the RAFs 
of galaxies more luminous than Mb = — 19 (roughly corre- 
sponding to Mk = —23.2) in groups of different richness 
(4 - 8, 9 - 20, 21 - 50, and > 50). We translate these 
richness classes into mass bins (log Af 2 oo = 13.78 — 13.90, 
13.90 - 14.36, 14.36 - 14.83, > 14.83), using the corre- 
lation between the number of galaxies and cluster mass 
presented in LMS04 (see Eq. IA1[) . Such a binning in mass 
is much finer than our simple low v.s. high mass sepa- 
ration, and reduces the size of cluster subsamples. Nev- 
ertheless, for log Pum = 23, we find the RAFs in these 
subsamples are 2.1, 3.1, 2.6, and 1.5%, respectively. Al- 
though the trend with respect to cluster mas s is in genera l 
consistent with that presented in Fig. 7 of IBestJ ( 20041 ). 
we do not find a strong suppression of RAF in the most 
massive clusters. Based on various mass binning and lim- 
iting radio and fC-band luminosities, we conclude that the 
RAF of the most massive clusters (e.g. log A/200 > 14.5) is 
slightly smaller than that in intermediate-to-massive clus- 
ters (e.g. log Af2oo = 14— 14.5), but the two are consistent 
within the uncertainties. In most of the cases, the RAF 
of the most massive clusters is greater than that of less 
massive clusters (e.g. logAT 2 oo < 14). 

A recent study of the optically-selected AGN population 
in clusters finds that the AGN fracti on is anti-correlated 
with the cluster velocity dispersion ((Popesso fc Bivianol 
20061). As the velocity dispersion is positively correlated 
with cluster mass, this finding appears to be an oppo- 
site trend to our results. However, as the radio and op- 
tical (emission-line) AGN activities may be independent 
of each other (see discussion in H5.4p . the distinct trends 
with cluster mass of the radio and optical AGNs probably 
point out to the different physical mechanisms that are 
responsible (e.g. mergers or interactions between member 
galaxies v.s. cooling instability). 



4.2.2. Comparison with the Radio Active Fraction in the 
Field Environment 

Combining the fraction / of galaxies that are radio- 
active as a function of their stellar mass M. with the stellar 
mass function (j>, one can estimate the RAF in the field en- 
vironment; that is, 

RAF - !KMWM)d\ogM 
f4>( M)d\ogM 
Based on Fig. 2 presented in Bes t et al.l ((2005) , we es- 
timate that, for hosts of logP > 23 AGNs, f(M) = 
1 ~ 28 ■ 79 ( M / M ) 2 ■ 43 . Using the stellar mass function de- 
termined bv lBell et al.l(|2003f ). and th e stellar mass-to- light 
ratio of 0.7 solar unit (if-band; see iLin et al.ll2003] ). we 
find that the RAF in the field is 1.3 ± 0.5% for M K < -24 
galaxies. For logP > 25, we approximate the result of 
iBest et all (I2005T) by f(M) = 10- i3 - 96 (M/M Q ) 3 - 63 , and 
find the fraction is 0.21 ± 0.08% for M K < -25 galaxies. 
These values are about 10 — 20% of the RAFs found in clus- 
ters (4.9 ±0.5% and 2.0 ±0.3%, respectively; see Table©. 
Together with the finding that the cluster RLF has higher 
amplitude than the (scaled) RLF in the field ( q4.1[) . these 
results suggest that cluster environment promotes AGN 
activity. 

4.3. AGN Duty Cycle 

The simplest application for the radio active fractions 
presented in the previous section is to estimate the duty 
cycle of radio-loud AGNs. If every galaxy has a central 
supermassive black hole (SMBH), the RAF is equivalent 
to the probability that the central SMBH is accreting for a 
galaxy of given luminosity range. In addition, if the AGN 
activity is intermittent and the RAF / remains about con- 
stant over the lifetime of the galaxies, the lifetime of a radio 
source is sim ply t r ~ ft„, w here t g is the typical age of 
a galaxy (e.g. ISchmidtj|1966l ). For a M* galaxy formed at 
z = 3, t g = 11.4 Gyr; then / « 0.05 gives t r w 5.7 x 10 8 yr. 
Based on bias and mass information from the spatial clus- 
tering of quasars, limits can be derived for the lifetime £q 
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of qu asars (e.g. iMartini fc Weinberdl200lt [1 aiman fc Huil 
20011); the recent determination from the 2d F QSO survey 
finds that 4 x 10 6 < t Q < 6 x 10 s yr. at z ~ 2 l|Croom et all 
2005), consistent with our estimate. At lower redshifts 
(0.05 < z < 0.95), the fraction of galaxies which exhibit 
AGN signature in their optical spectra may be as high as 
40% suggesting a duty cycle of w 2 x 10 s yr. (|Miller et alj 
2003), which is within the uncertainties from our value. 

The higher RAF of BCGs strongly suggests that their 
duty cycle is longer than the other galaxies, which may 
be closely related to their location in clusters and the am- 
ple supply of gas due to radiative cooling of the ICM. We 
explore the implication of our results on the ICM energy 
injection from the radio-loud AGNs in £j6j 

4.4. Comparison with Results of \Ledlow & Oweri 

In Sj4T]we compare the number densities obtained from 
the cluster and field RLFs, and find that the clusters con- 
tain more radio-loud AGNs than the field does, even after 
the difference in the overdensities has been accounted for. 
On the other hand, L096 cast their results in terms of frac- 
tional luminosity function (meaning that the RLF is fur- 
ther normalized to the total number of galaxies expected 
in each radio luminosity bin), which can be compared to 
the RAFs presented in tj4.21 

L096 consider both the univariate and bivariate LFs (in 
our notation, ip(P,M R < -20.5) and ip(P,M! < Mr < 
M-i), respectively). Summing up the univariate LF (pre- 
sented in their Fig. 6), they find that RAF=14% (for 
log Piim = 22) for early type galaxies with Mr < —20.5, 
which can be compared to our value of 12.5 ± 3.0% (for 
MK.Um — —23.1, where R — K = 2.6 for early type galax- 
ies is used) , although we do not distinguish between early 
and late type galaxies. 

Next we compare our data with their bivariate LF for 
galaxies with i?-band magnitudes within —23.5 < Mr < 
-22.7 (roughly corresponding to -26.1 < M K < -25.3). 
At logP = 23.8, 24.2, 24.6, and 25, the RAFs (in percent- 
age) based on their Table 2 and our data (in parantheses) 
are: 2.2 ± 0.9 (4.0 ± 0.9), 4.0 ± 1.2 (4.1 ± 0.8), 5.8 ± 1.5 
(1.7 ±0.5), and 1.8 ±0.8 (1.5 ±0.5). 

There are several differences that make the comparison 
between the two analyses challenging, including: (1) the 
distinction of galaxy morphological types, (2) the portion 
of the cluster surveyed, (3) the deprojection of the RAF 
from a cylindrical to a spherical volume. Given the relative 
abundance of the early and late type galaxies in clusters, 
and the fact that radio-loud AGNs are mainly hosted by 
early types, exclusion of the late type galaxies should not 
be a problem. However, the second and third factors do 
affect the calculation of the RAF. The RAF will be radius 
dependent because of the rather different radial distribu- 
tion of radio sources and galaxies in clusters (see Fig. 01 . 
Let fs denote the RAF measured at a radius correspond- 
ing to an overdensity 5, one expects that fs 1 > f$ 2 for 
Si > 62- L096 calculate the RAF for galaxies within 0.3 
Abell radius (about 0.6 Mpc). This fixed metric radius 
corresponds to different overdensities for clusters of differ- 
ent masses, and the RAF from L096 should be regarded as 
an averaged value; without specific knowledge of the mass 
distribution of their cluster sample, it is not possible for a 
direct comparison. However, unless the majority of their 



clusters are low mass systems, we would expect their RAF 
to be higher than ours, given that 0.6 Mpc is within r 50 o 
for a moderately massive cluster (e.g. M 200 ~ 5xlO 14 M ). 

Furthermore, we expect their RAFs to be increased by 
about 30%. This is because they do not explicitly take 
into account the projection of galaxies along the line of 
sight. Knowing the radial profiles for galaxies and radio 
sources enables one to calculate the correction factor to 
transform the RAF measured within a cylindrical volume 
into one that is within a sphere of the same radius. 

LQ96 compare their RLF with that obtained by 
ISadler et all 0989) and find no significant difference be- 
tween the RLFs measured in the field and in clusters. How- 
ever, because of the reasons mentioned above, it may not 
be a fair comparison; if their RAFs were to be increased 
(i.e. due to deprojection), then the cluster RAF would 
be higher than the field value. Furthermore, as the two 
measurements are made at different frequencies (1.4 and 
5 GHz), the effect of spread in the spectral index should 
be taken into account f £|2.4[) , which makes it harder to 
precisely compare the results. 

To summarize, a more direct comparison with the re- 
sults of L096 is made possible with our XC sources. Al- 
though agreement is found for both univariate and bivari- 
ate LFs between the two groups, we point out differences 
in the analyses (the fraction of the cluster region surveyed 
and the deprojection of the RAF) that may make their 
RAFs systematically higher than ours. A possible expla- 
nation of the agreement is that the cluster sample used 
by L096 is on average less massive than ours, so that the 
effects from the two factors mentioned above (portion of 
the clusters surveyed and the deprojection) roughly cancel 
out. 

5. ACTIVE RADIO SOURCES IN LOW-MASS CLUSTER 
GALAXIES 

The higher amplitude of the RLF in the cluster versus 
the field is an exciting result that deserves further con- 
sideration. We first discuss the possibility that the lower 
amplitude in the field RLF is due to a different selection 
technique, then examine the robustness of our determina- 
tion of the cluster RLF. 

To determine a field RLF, a usual procedure is to first 
define a sample of galaxies limited by, say, an apparent 
magnitude mu m , then conduct the radio survey for these 
galaxies, and follow up with optical s pectroscopy for the 
sources that are cross-matched (e.g. ISadler et al.1 [l989; 
ICondon et alll2002t ISadler et al1l2002fl . We note that in 
doing so, the radio galaxy sample is only complete to the 
absolute (optical) magnitude that the highest-redshift ob- 
ject in the sample corresponds to. Such a survey will there- 
fore miss the radio sources in low luminosity galaxies, and 
the LF thus constructed will underestimate the true am- 
plitude and shape of the RLF. In some sense, the field 
RLF is built from objects like our XC sources. The com- 
parison, therefore, would be more meaningful between the 
field RLF that is based on sources selected with a fixed 
OIR absolute magnitude limit and the RLF from our XC 
sample. 

Of course, if powerful radio sources (e.g. logP > 24) 
only reside in luminous galaxies, then our concern for the 
field RLF will be a non-issue. We will come back to this 
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point below. 

We show in Fig. [8] as the hollow symbols the cluster 
RLF based on the XC sources that have K-b&nd mag- 
nitude M K < -24 (~ M K *] e.g. LMS04) in 202 clusters. 
The highest redshift of clusters that we can include so that 
the galaxies are complete to Mk = —24 is z = 0.0693. The 
background contribution is accounted for as follows: given 
a limiting i-T-band absolute magnitude, for each cluster 
we calculate the corresponding apparent magnitude Kf, 
and use the (radio) log N— log S from the XC sources with 
K<K f . 

The RLF for all NVSS sources from all our clusters is 
shown as the solid points in Fig. [5] The difference in ampli- 
tudes between the NVSS-only and the XC RLFs encodes 
the fraction of all radio sources that are hosted by lumi- 
nous, presumably massive, galaxies, as the i-T-band light 
traces stellar mass well. For example, down to P = 10 23 
WHz -1 , 55 ± 14% of all radio sources are within galaxies 
more luminous than Mk — —24. Repeating the exercise, 
we find that 43 ± 11% (52 ± 15%) are in galaxies with 
Mk < —25 (—23). From these comparisons we find that 
the fraction of cluster radio sources associated with lumi- 
nous galaxies is ~ 1/2. 
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Fig. 8. — Comparison of RLFs constructed from the NVSS-only 
(solid symbols) and XC (hollow symbols) sources. The crosses show 
the RLF from a smaller sample of those clusters that have redshift 
information for their member radio galaxies (based on data from 
IMiller fc O"wer]|200ll) . 

Earlier surveys conclude that powerful radio sources are 
almo st exclusively found in luminous galaxies (e.g. Mr < 
-21: ISadler et al.lll989t lOwen fc Whitdll99H LQ96). In 
particular, L096 notice that all galaxies fainter than 
Mr = — 21 that are identified with radio sources stronger 
than the flux limit of their survey turn out to be back- 
ground objects. First, we note that Mr, = —21 roughly 
corresponds to Mk — —23.6 (for early-type galaxies), only 
~ 0.5 mag fainter than the characteristic magnitude of the 
luminosity function. The i-T-band data of our XC sources 
(more powerful than e.g. logP = 24) reach deeper than 
that limit (e.g. Mk > —23, although the sample size is 
smaller; see Fig [B] above) . It is not clear to us either what 
(optical) magnitude limit the L096 survey reaches, or the 



number of faint galaxies with redshift measurements that 
are determined as background. 

As another check o f our RLFs, we use the data from 
IMiller fc Owed (|200lD , who have conducted an extensive 
redshift survey over 18 nearby (z < 0.033) clusters and 
provided the redshifts for a sample of 326 radio galax- 
ies. Eight of these clusters are in our sample (so that 
we can estimate their masses), and we cross-match their 
source list with the 2MASS catalog. The RLF from this 
dataset is shown in Fig. [5] as crosses (for 20 galaxies with 
Mk > —24). It shows good agreement with our XC RLF 
within the (rather large) uncertainties. 

5.1. Robustness of Our Treatment of the Radio 
Luminosity Function 

The analysis presented above indicates that about 50% 
the cluster radio galaxies in our sample are associated with 
galaxies less luminous than Mk = — 23. Because this is a 
new result, we discuss some possibilities that may explain 
the differences between the RLFs of the NVSS-only and 
XC sources. These include: (1) differential breakdown of 
the statistical background subtraction scheme, presumably 
due to the differences in the dominance of the background 
(the total-to-background number ratio is > 7.6 and 1.3 for 
the XC and NVSS-only cases, respectively), (2) presence 
of radio sources that are resolved into several components, 
(3) missing matches due to astrometric problems and (4) 
gravitational lensing of background sources. 

For the first possibility, we have examined (a) overall 
shape and amplitude of the logA^-logS; (b) variation in 
the background as a function of declination 8 or Galactic 
latitude b (see ij2.2p : and (c) the effects of sidelobes. We 
find that our log iV-log S has to be increased in amplitude 
by ^ 10% to account for the differences in RLFs, which 
is unlikely, for it is in good agreement with that derived 
from the entire NVSS survey region (6 > —40 deg). In 
addition, the mean surface density thus de rived is in ex- 
cellent agr eement with the val ue obtained bv lBlake fc Wall 
(|2002f) and lBlake et alj (|2004l ). which is derived after care- 
fully masking out from the NVSS catalog the Galactic 
plane, nearby bright extended sources, and sidelobes of 
very strong sources. Secondly, we do not find any clear 
changes in amplitude for both NVSS-only and XC RLFs 
as a function of either b or S. Thirdly, to minimize the ef- 
fects due to sidelobes that are not cleaned from the NVSS 
catalog, we construct RLFs from clusters that do not have 
sources with flux density larger than 1 Jy. The RLFs are 
statistically indistinguishable. 

Finally, the RLFs are not sensitive to our choice of the 
completeness limit. For example, using 15 or 20 mJy as 
the limit does not change either the shape or amplitude of 
the RLF for bins at log P > 23.5. 

Although there has been an estimate of the fraction of 
NVSS sources that are c omponents of the same source 
(~ 7%, iBlake et al.ll2004fl . their effect on the RLF is not 
trivial to evaluate. To examine if the resolution of the 
NVSS affects the resulting RLF, we construct RLFs both 
for radio-source only and for radio-2MASS cross-matched 
sources, using the data from the FIRS T (Faint Images o f 
the Radio Sky at Twenty-centimeters, iBecker et al.l fl995) 
survey, which has a resolution of 5". Because of the sky 
coverage of the survey, only 242 clusters are available. We 



Radio Galaxies in Clusters 



15 



use only sources with S > 10 mJy. The FIRST-only RLF 
is shown in Fig. [5] as hollow squares, and overall it is in 
agreement with the NVSS-only RLF. We also find that 
the XC RLF based on the FIRST survey is consistent 
with that determined using the NVSS-2MASS XC sources. 
Integrating down to logP = 23, the FIRST-2MASS XC 
sources has the space density that is 0.98 ± 0.23 of that 
from the NVSS XC sources. Therefore, there is no evi- 
dence that the resolution of the NVSS survey is causing 
problems. 

Astrometry should not be a problem, as 2MASS ex- 
tended sources have accuracy good to < 0.5", and NVSS 
sources with S > 10 mJy are accurate to 1 — 2", which are 
much smaller than our maximal separation d max = 20" 
for cross-matching the two cata l ogs. 

In the case of lensing, [Cooravl (|1999f ) estimates the num- 
ber of expected lcnsing events due to clusters more massive 
than ~ 1.3 x 10 15 M Q over the sky, and find it to be ~ 1500 
for sources with S > 10 /iJy at 1.4 GHz. The lensing prob- 
ability p depends on both th e number dens ity n and mass 
M of clusters (p ~ nM 4//3 ; ICooravlll999f h Considering 
the shape of the cluster mass function, and our adopted 
flux limit (10 mJy), lensing events due to our c luster s 
may not be larger t han the prediction of ICooravl jl999). 
iPhillips et all (|200l[ ) conduct a survey over 4150 deg as 
an attempt to detect lensed quasars by clusters of masses 
> 10 13 - 10 14 M Q . None of their 38 candidates are found 
to be actual lensing systems. Therefore, we conclude that 
the frequency of lensing event s is quite small (ICoorav et alj 
119981: lAndernach et alJll998t see also iBovce et al.ll2006ft . 

None of these alternatives seem to be able to offer a 
satisfactory explanation. Below we present some evidence 
that supports a population of AGN in low mass galaxies. 
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Fig. 9. — Total surface density profile (cluster plus background) 
for the NVSS-XC sources (the NVSS sources without a 2MASS 
counterpart within d m ax chosen for each cluster, see £|2.2I hollow 
squares). The profile for all the NVSS sources is shown as solid 
points. Both are from sources with log P > 23. The two profiles are 
arbitrarily normalized to match at the innermost bin. The fact that 
NVSS— XC sources are clustered, together with the small expected 
effects from gravitational lensing amplification, indicate that they 
are very unlikely to be background sources. 



5.2. Evidence Supporting the Existence of Such 
Population 

We can estimate the RAF expected for the low mass 
cluster galaxies so that the difference between the NVSS- 
only and XC RLFs can be accounted for. The basic idea is 
to calculate, for galaxies whose radio luminosity is above 
a given limit, the ratio between the numbers of galaxies 
more and less luminous than a certain AT-band magnitude 
(which we choose to be Mk = —23). For the KLF with the 
faint-end power-law slope a — —0.9, for galaxies within the 
magnitude range from —23 to Mu m , where Mu m = —19, 
-18, and -17, the RAFs are 1.84%, 1.51%, and 1.30%, re- 
spectively. Therefore the active fraction of these low-mass 
cluster galaxies does not need to be high given reasonable 
forms of the KLF. 

Further support of our conclusion comes from the spatial 
distribution, as well as the RLF of these sources. We first 
examine the surface density profile for the NVSS sources 
that do not have a 2MASS counterpart (hereafter referred 
to as the NVSS-XC sources), and find that these sources 
are clustered (Fig. [3]). If they were all background objects 
we would not expect them to be clustered around the po- 
sitions of known galaxy clusters. 

Secondly, we construct the RLF for the NVSS-XC 
sources. For this purpose we use the log A-logS 1 based on 
NVSS-XC sources from a region of 60 deg radius centered 
at the NGP. We follow the same procedure as outlined 
in H4.ll and the resulting RLF from all clusters is shown 
in Fig. [TUJ Integrating the two RLFs down to logP = 23 
(24), we find that the number density of NVSS-XC sources 
is 67 ± 18% (55 ± 11%) of the NVSS-only sources, which 
is similar to the estimate based on the NVSS-only v.s. XC 
RLF comparison (jjSJ). 

Is there evidence from OIR observations for the exis- 
tence of low mass cluster galaxies that host radio-loud 
AGNs? We have obtained AT-band imaging data sufficient 
to probe to Mk ~ —19 (about five magnitudes fainter 
than Mk*) for a sample of 10 clusters at z < 0.2 (Lin 
et al. 2006, in preparation). Due to the field-of-view of 
the observations (20' x 20', conducted with FLAMINGOS 
imager at KPNO 2.1m), we focus on the region enclosed 
by the radius T2ooo (within which the mean overdensity is 
2000p c ) for this sample of clusters. There are 21 sources 
in the NVSS catalog that are within r 2 ooo in these clusters 
and are stronger than 10 mJy; we expect 4.44 background 
sources based on the mean surface density of the NVSS. 

Because all 10 clusters lie within the region covered by 
SDSS data release 4 (DR4), we seek optical and A"-band 
counterparts of the NVSS sources. We find matches in 
both the A-band and optical catalogs for all (21) but 
one source (down to our completeness limit K ~ 19.3). 
In Fig. [11] we show the g — r v.s. r color-magnitude di- 
agrams for both the radio sources (hollow symbols) and 
all galaxies (small points) projected within r^ooo- In all 
cases the red sequence of the cluster galaxies is apparent. 
As a guidance for cluster membership we gather redshift 
information from both SDSS and the NASA/IPAC Ex- 
tragalactic Database (NED) and show the galaxies whose 
redshift is within 0.001 from the cluster redshift as the 
large solid points. It is clear that most of the optically 
bright (r < 19) radio sources are cluster members. Using 
the red sequence as a rough guide for the cluster member- 
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ship assignment, and taking into account the uncertainties 
in the photometry for the faint sources, it is likely that 
w 1/2 of the optically faint (r > 20) radio sources are 
members. Therefore, we conclude that the background es- 
timates from both the NVSS log iV-log S and the optical 
color of the XC sources give roughly consistent results (the 
former method gives 16 — 17 members, and the latter yields 
an estimate of 14 — 15 members). 




22 24 26 

log[P 1 .4GH Z /( w Hz" 1 )] 

Fig. 10. — RLF for NVSS-XC sources (hollow squares). As a 
comparison, the RLF for all NVSS sources is shown as the solid 
points. The RLF of NVSS— XC sources appears to be similar to 
that of NVSS— only sources. 

We examine the radio luminosity for candidate cluster 
member galaxies of r ~ 20 — 21 in four clusters (A85, 
A2255, A2055, A1914), and find that P = 4x 10 23 - 2 x 10 24 
WHz -1 . These luminosities are slightly smaller than that 
where the discrepancy between the NVSS-only and XC 
RLFs is largest (e.g. logP = 24 - 25, see Fig. E)- We 
expect, however, if a similar analysis is carried out for a 
large, nearby cluster sample (e.g. 200 clusters), we will find 
low-mass radio galaxies whose radio luminosity is greater 
than those found here. 

5.3. Seeking Spectroscopic Confirmation 

Conclusive support for the existence of the active low- 
mass cluster galaxies would come from spectroscopy. From 
Fig. [Tl]it is clear that the majority of such sources, typi- 
cally having r < 20, would lie below the limits of current 
large galaxy redshift surveys (e.g. SDSS and 2dFGRS). 
Nevertheless, we combine the three large public datasets 
(NVSS, 2MASS, SDSS DR3) to seek further evidence. In 
Fig. [12] we show the distribution of radio sources in the 
Mk-P space, based on a sample with measured redshifts 
from the SDSS survey. More specifically, we cross-correlate 
our XC sample (not limited to those within cluster fields) 
with the main galaxy sample in the SDSS DR3. There are 
1318 galaxies matched that are above the completeness 
limits of both 2MASS and NVSS. 

We caution that the majority of these galaxies would 
not be cluster galaxies. More importantly, we have not 
assessed the statistical completeness of the sample, as well 
as taking into account the variation of the volume probed 



at different locations on the Mk-P plane. The plot is to be 
used as merely a suggestion for the form of the BLF; the 
concentration of data points on the figure should not be 
taken to be proportional to the actual BLF. Despite these 
caveats, the figure reveals a couple of interesting proper- 
ties of galaxies in the local (z < 0.2) Universe. For massive 
galaxies (e.g. Mk < —24), at a given Mk, their radio lu- 
minosities span a very wide range. This is to be expected, 
for at the moment we observe, the central SMBHs of galax- 
ies must be at various stages of their duty cycle. Secondly, 
the maximum radio luminosity the galaxies can have seems 
to correlate with th eir stellar mas s, which is in agreement 
with the findings of I Owed (|l993f ) and L096. This result 
can be expected from the correlati on between the masses o f 
the bulge and the central SMBH (|Magorrian et al.lll998l) . 
and the classical idea of maximum luminosity given the 
mass of the source (the Eddington limit; e.g. Ide Young] 
120021) . 
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Fig. 12. — Correlation between 1.4 GHz and ii"-band luminosities 
for (predominantly field) galaxies at z < 0.2. We cross correlate 
our XC sample with the main galaxy sample from the SDSS DR3; 
in total there are 1318 matches that are above the completeness 
limits of 2MASS and NVSS. The density is shown in logarithmic 
scale, where the units are arbitrary. We caution this sample is for 
illustrating the shape of the bivariate LF only, and is not meant for 
rigorous statistical analysis. 

If we restrict the galaxies to lie within clusters, the 
resulting distribution is similar, but with far fewer data 
points. This simply indicates that the SDSS is not deep 
enough (and is not designed) to probe such a faint popula- 
tion. We next search for evidence from two large surveys 
targeting cluster galaxies: the E SQ Nearby Abell Cluster 
Survey (ENACS. lKatgert et al1ll996l ) a nd the Cluster an d 
Infall Region Nearby Survey fCAIRNS. lR.ines et al.ll2003t) . 
By cross correlating the NVSS catalog with these surveys, 
we do not find any faint cluster galaxies (e.g. Mk > —23) 
that are associated with strong NVSS sources. We note 
that, however, for our purpose, the depth of these sur- 
veys is also not sufficient. For example, the completeness 

of the CAIRNS is well below 50% at about M K 21. 

Moreover, our estimate is that only 1 to 2% of these faint 
galaxies are radio luminous, and so we require large num- 
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Fig. 11. — Color-Magnitude diagrams for a sample of clusters covered by SDSS, for which we have obtained deep iC-band imaging data. All 
SDSS extended sources projected within r2000 from cluster centers are shown as small points. For a limited number of galaxies whose redshift 
is available (from SDSS and NED) and is within 0.001 from the cluster redshift, they are displayed as large solid points. The NVSS sources 
with iC-band and SDSS counterparts are shown as hollow symbols with errorbars. Among these sources, those that are faint (r > 20) and lie 
along the cluster red-sequence are candidates of the NVSS-XC sources whose optical counterpart is too faint to be detected by 2MASS. 



bers of faint galaxy spectra to be assured of a single con- 
firmed AGN. Deep spectroscopy extending to Mk ~ — 19 
is needed in a dozen or so clusters to examine the low mass 
cluster radio galaxy population. 

5.4. Discussion 

A recent comprehensive study of the properties of radio 
galaxies in the local Universe has shown that the fraction 
of galaxies that host radio-loud AGNs is a strong func- 
tion of both the BH and stellar masses (|Best et all 120051 
hereafter B05). To be more specific, they find that RAF 
cx M s 2 4 and RAF oc M BH , where M B h and M sta r are su- 
permassive black hole and stellar masses, respectively. The 
study also shows that the radio and (optical) emission-line 
AGN activities are independent of each other. B05 con- 
clude that the radio selection picks up the most massive 
SMBHs that are believed to be largely dormant at present. 

Can this strong dependence on stellar mass of the nu- 
clear radio activity be reconciled with our proposed pop- 
ulation of radio-active low-mass cluster galaxies? A possi- 
bility lies in the confining pressure that the ICM imposes 
on cluster galaxies. For an elliptical galaxy surrounded 
by hot gas that is isothermal with temperature T and has 
the X-ray luminosity Lx , we follow the simple arguments 
outlined in B05 and estimate the rate that gas cools out to 
be Mcooi oc Lx/T. Based on the observed scaling relation 
between Lx and optical luminosity of galaxies, the Faber- 
Jackson relation, and that between Mbh and the velocity 
dispersion a, B05 argue that 

M coo i oc cr 6 oc Mbh- 



We note that because the galaxy sample B05 use is drawn 
from the main galaxy sample of the SDSS, they do not 
p robe to the low mass re g ime in cluster galaxies. 

iFinoguenov fc Miniatil (|2004h examine the effects of 
cluster environment have on the Lx-L opt a 2 relation 
(where L opt is the optical luminosity of galaxies) and find 
that it flattens in clusters with respect to the field. The 
inferred relation in cluster is Lx oc J lll / 4 j which leads to 

Mcooi oc cr 7/2 

(using T oc cr 2 for an isothermal gas). Therefore, if the 
Mbh — & scaling remains the same, the cooling rate be- 
comes only weakly dependent on stellar or black hole mass. 
Although we do not consider several other factors such as 
the radiative accretion efficiency onto the central SMBH 
and the abundance of the gas that feeds the SMBH, the 
change in the gas cooling rate suggests that a simple ex- 
trapolation of the RAF from properties of high mass galax- 
ies is unlikely to hold. Furthermore, we note that the RAF 
is a strong function of the radial distance from the clus- 
ter center. There is clearly something special about being 
deeper within the cluster potential well- at least for the 
high mass galaxies, and so it is reasonable that this physics 
is also at work for the low mass galaxies. More rigorous 
investigations, both theoretical and observational, are def- 
initely needed to further study this population. 

6. ICM ENERGY INJECTION FROM AGNS 

We start by summarizing several observational results 
about the BCGs and the central region of clusters. (1) The 
spatial distribution of radio-loud AGNs is much more cen- 
trally concentrated than both the ICM and cluster galaxies 
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(©. (2) About 80% of the BCGs are foun d within 0.1r 200 
from the cluster center (|Lin fc Mohrll2004f) . (3) BCGs are 
much more likely to host AGNs of a given radio luminos- 
ity compared to less massive galaxies ( ij4.2p . (4) BCGs 
have 2 — 3 times higher likelihood to host AGNs compared 
to the non-BCG galaxies of comparable luminosity. (5) 
The BCGs in high mass clusters show a ~ 3 times higher 
probability of hosting an AGN compared to those in lower 
mass clusters ( 34.20 , 

It is clear that the BCGs spend more time in the ac- 
tive phase than the other cluster galaxies do. This also 
implies this elevated AGN activity is linked to the en- 
vironment of BCGs; apparently the cl uster core is a fa- 
vorab le place to turn on an AGN (e.g. lBurndll99d: lEilekl 
2004). This is most likely due to the repeated merg- 
ers central galaxies exp erience over the a ssembly his- 
tory of the cluster (e.g. iLin fc Mohrl [2004). and/or in- 
teractions with the dense ICM surrounding them. All 
the observations mentioned above favor the idea that 
radio activities of central galaxies may be a plausible 
source to halt the radiative cooling of the ICM toward 



the center (e.g. Binnev fc Tabor 19951: iDavid et alj 
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With the recent finding that the radio luminosity at 1.4 
GHz may be a rough indicator of the mechanical luminos- 
ity L mec h the AGN outflow ("bubble s" or "cavities" see n 
in X— ray images) has on the ICM (|Birzan et~all[200l . 
we can use the AGN RLF to estimate the amount of en- 
ergy injected into the ICM by the cluster AGN popula- 
tion, and comp are it to the t h erma l energy content of the 
ICM. Following iBirzan et all (|2004f) . who find that L mech 
is of order 100^P„ where v = 1.4 GHz, we assume that 
Lmech — AQQi]vP Vl where 77 is an efficiency coefficient, and 
the factor of 400 comes from our further assumption that 
the gas co ntained in the AGN outflows is relativistic and 
adiabatic (jBirzan et al.ll2004T ). With this expression, we 
can calculate the total energy output from AGN Eagn 
by integrating the mechanical luminosity over time from 
z = to z = 1. A 7 = 2.5 power-law evolution in AGN 
population is also assumed (see §T]). 

The ICM thermal energy is estimated as E t h — 
{3/2)kTxfiCMM 200 /^m p , where f W M = 0.12 is the ICM 
mass fraction (assumed to be constant with respect to clus- 
ter mass), = 0.59 is the mean molecular weight, and m p 
is the proton mass. For a M200 = 10 15 M Q cluster, we 
find that within the virial radius, Eagn / 'Eth = 0.0047?7; 
for a low mass cluster (10 14 M Q ), the ratio is 0.02047?. Be- 
cause of the different mass dependences of Eagn and Eth, 
the AGN heating is more relevant in low mass systems 
than in massive ones. We can also express the amount 
of AGN heating in terms of the energy per ICM particle, 
which would be independent of cluster mass. We find that 
Eagn corresponds to 0.05977 keV per particle. 

These numbers suggest that AGN mechanical heating 
is not a significant contributor to the thermal energy of 
the ICM within the entire cluster virial region. How- 
ever, while the ICM is rather extended, the AGN dis- 
tribution is very centrally peaked; this implies that the 
AGN would be more energetically important toward the 
central parts of clusters. We use the measurements of 
the gas density profile for a sample of 45 nearby clusters 



(Mo hr et al.| [i999) to determine the ICM thermal energy 
within 10% of the virial radius. We find that within this 
radius, E AGN /E th = O.OII77 and 0.046?7 for 10 15 M Q and 
10 14 M Q clusters, respectively. Cast differently, the AGN 
energy is equivalent to 0.137? keV per particle. It is ap- 
parent that, as long as the efficiency coefficient is not too 
small (e.g. 77 > 0.1), the AGNs are an imp ortant source 
of hea ting of the ICM. Based on the data of IBirzan et alj 
(2004]), 77 may range from 1/30 to 20. Another complica- 
tion is that, 77 may be strongly time-varying, which makes 
the amount of heating extremely difficult to estimate. Our 
crude calculation is only meant to be an order of magni- 
tude estimate of the role AGNs can play in the thermody- 
namic history of clusters. 

7. IMPLICATIONS FOR SZE SURVEYS 

Our results can be used to estimate the contamination 
of the cluster SZE signal due to radio point sources. We 
address this issue in two steps: we first transform the ob- 
served 1.4 GHz AGN RLF into the frequency at which an 
SZE experiment operates; in the second step we use a sim- 
ple Monte Carlo scheme to estimate the fraction of clusters 
whose AGN flux is a significant fraction of the SZE flux. 

Given the frequency v (in GHz) of an SZE experiment, 
we transform our observed 1.4 GHz RLF to that frequency, 
via 



dn 



d\ogP v 



(log P.) 



dn ( . 1.4 \ ,. , , 

— — a log V log P v f(a)da 

alogPi.4 V v J 

(4) 



where P v and P\a are the luminosities at v and 1.4 
GHz, respectively, which are related to each other by 
Pi a = (lA/v) a P v ; /(a) is the probability distribution 
of the spectral indices from cluster radio sources, taken as 
the SID presented in 32.41 (see Fig. [3]). (We note in pass- 
ing that ILin et al.ll2002l use an expression very similar to 
Eqn. 2] to extrapolate the LF to higher frequencies.) To 
account for the possible steepening of the spectral index at 
high frequencies (e.g. v > 100 GHz), the LFs at frequen- 
cies greater than 90 GHz arc obtained by extrapolating 
the 90 GHz RLF with a modified SID f{a) = f(a + 0.5). 

In Fig. [TBI we show the estimated AGN RLFs at 15, 30, 
90, and 150 GHz. The most apparent feature in these ex- 
trapolated RLFs is the rise of the bright end, which is due 
to the non-negligible population that has a rising spec- 
trum (a > 0; see iJ2.4|) . and the steepness of the RLF. As 
the bright end will definitely contribute significantly to the 
fluxes from radio sources, this points out the importance 
of an accurate model for the SID of the radio sources over 
the frequency range of interest. 

Our technique involves an extrapolation by as much as 
a factor of 100 in frequency, and so this simple approach of 
using the 1.4 GHz to 4.85 GHz spectral behavior together 
with a break at 100 GHz should be viewed as a crude 
extrapolation- especially at the highest frequencies. Ex- 
isting observations tend to suggest a flattening of the spec- 
tral index for samples selected at higher frequencies. For 
example, a VLA survey finds that for sources (mainly flat- 
spectrum, core-dominated QSOs) with S > 0.1 Jy, d(8, 90) 
can be described by a Gaussian wh ose mean and dispersio n 
are -0.37 and 0.34, respectively (|Holdawav et al.l ll994D . 
At fainter flux limit (~ 20 mJy), results from the CBI 
experiment provide an estimate that d(1.4,31) = —0.45, 
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with maximu m and minimum ind exes being 0.5 and —1.32, 
respectively ( Mason et al.l l2003f) . An extensive follow- 
up from 1.4 to 43 GHz of the 9C survey finds that all 
but two sources in t heir sample of 176 show a falling 
or concave spectrum (jBolton et al.ll2004t) ; about 20% of 
their sources have spectra that peak at or above 5 GHz. 
Wit h sensitive (0.4 mJy) 1cm observations toward 56 clus- 
ters, ICooravet all (|l998T l find a(1.4, 28.5) = -0.77 ± 0.06 
from 54 sources. A recent follow u p study finds tha t 
5(1.4,28.5) = -0.70 from 87 sources (jCoble et al.l 12006ft . 
Finally, the WMAP satellite has conducted an all-sky sur- 
vey of bright radio sources (> 1 Jy) at 41 — 94 GHz, and 
found that the mea n index a w with a scatter ~ 0.3 
(jBennett et al.|[2~003l ). It is clear that reliable estimates of 
the cluster radio galaxy luminosity function at high fre- 
quencies will have to await new observations. 
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Fig. 13. — The AGN luminosity functions (within r20o) at several 
frequencies. We transform the 1.4 GHz cluster AGN RLF (solid line) 
to 15, 30, 90, and 150 GHz (dotted, short-dashed, long-dashed, and 
dot-dashed lines, respectively), using Eqn. [4] Also shown as points 
is an estimate of the 30 GHz c luster radi o source LF, based on the 
data presented in lCoble et al.l (120061 ') and lCooray et al.l (119981). Th e 
background is estimated using the log jV-log S of lKnox et al.I f2004'). 
The agreement between the points and the short-dashed line implies 
our method of extrapolating the RLF works well. 

The data points with uncertainties in Fig. [T3] show an 
estimate of the cluster radio source LF within r2oo at 
30 GHz, ba s ed on 68 c lusters. The d ata a re taken from 
ICoble et al.l (|2006f) and lCooray et all (|1998fh Taking into 
account the beam response as a function of radius, we es- 
timate ithj3_bjJ£kground using the log iV-log S at 30 GHz 
(from [Knox et al.l I2004T ) . Despite a large degree of scat- 
ter, the match between the amplitudes of the data points 
and the short-dashed curve suggests that our extrapolated 
RLF at 30 GHz is in reasonable agreement with the avail- 
able data. 

For ease of comparison with other work, we have fit our 
extrapolated RLFs and provide the fitting parameters in 
Table [6] The same functional form as Eqn. [3] is used. 

Given the mass and redshift of a cluster, and a model 
of AGN population evolution (see below) , the RLF can be 
integrated to give the total AGN flux. We note, however, 
the flux thus obtained is usually dominated by the rare, 



very luminous sources, and therefore the effects of contam- 
ination from AGNs would be overestimated. To obtain a 
more realistic estimate, we employ a simple Monte Carlo 
simulation scheme. First, we integrate the RLF (from 10 20 
to 10 28 WHz -1 ) to obtain the expected total number of 
AGNs, (N). Using (A'') as the mean, we generate a Poisson 
random number N p , to account for cluster-to-cluster vari- 
ations of the total AGN number. We assign fluxes to the 
Np sources, using the RLF as the probability distribution. 
The sum of the fluxes from these sources gives the total 
AGN flux (Sagn) of the cluster, which can be compared 
with the cluster SZE flux within the virial radius (Ssze)- 
We use the expre ssion of the SZE flux-clus ter mass re- 
lation provided in iMaiumdar fc Mohrl (|2004l ) . Repeating 
this process 10 5 times, we estimate the AGN contamina- 
tion fraction by calculating the fraction of clusters whose 
AGN flux is a significant proportion of their SZE flux. Let 
us denote Sagn > s\Ssze\- We regard a cluster as "lost" 
from a survey when s = 1 , and correspondingly define the 
"lost cluster fraction" (LCF) as the proportional of clus- 
ters whose AGN flux is equal to or greater than the SZE 
flux. 
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Fig. 14. — The cumulative fraction of clusters whose AGN flux is 
greater or equal to s\S$ze\ as a function of s. We consider the case 
for clusters of mass 2 X 10 14 Mq, at z = 0.6, measured at 150 GHz. 
The dashed line is the estimated fraction when the SID from all the 
sources for which we have flux measurements at both 1.4 and 4.85 
GHz are used to extrapolate the RLF, and the dot-dashed curve is 
obtained when the SID of the BCGs is used. The lost cluster frac- 
tion is the value of the curve at s = 1 (about 12% and 2% for the 
two cases). While the lower curve represents an extreme case, the 
upper curve should be more representative of all the sources that 
contribute to AGN contaminations. 

For our modeling of the AGN population evolution, we 
consider two cases: a no-evolution model (i.e. the AGN 
RLF is independent of z), and a power-law evolution model 
in which the number of AGN in a cluster of mass M scales 
as (1 -I- z) 7 , but the shape of the RLF remains the same. 
The true evolution towards high redshifts is still under de- 
bate. For example, an earlier study suggests that the AGN 
RLF in interme diate redshift clust ers is si milar to that in 
nearby clusters (|Stocke et al.lll999t see also lPerlman et ail 
2004| ICoble et al.ll2006l ). However, we note that the stud- 
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Table 6 
Luminosity Function Fits 



frequency 


y 


b 


X 


w 


1.4 


37.97 


2.40 


25.80 


0.78 


15 


40.31 


4.27 


27.10 


0.76 


30 


40.60 


3.91 


27.89 


0.81 


90 


41.53 


3.61 


29.12 


0.84 


150 


40.96 


3.34 


28.88 


0.85 



Note.- See Eqn. [3] for the functional form used in the fit. 



ies of IStocke et al.l and ICoble et al.l focus on radio-loud 
AGNs within a fixed metric radius over the redshift range 
of their cluster samples (0.3 < z < 0.8 and 0.1 < z < 1.0), 
and therefore the AGN populati on may not be prope rly 
compared (see also discussions in lBranchesi et al]|2006l for 
other possible explanations) . Studies that favor a positive 
evolution also exist. A deep radio survey over the cluster 
MS1054-03 (z = 0.83) finds that the radio-loud AGN frac- 
tion is ~ 4 times higher tha n that in z < 0.1 clusters found 
by L096 ({Best et alj|2002f) . Another determination of the 
RLF from a sample of 18 clusters at 0.3 < z < 0.8 also 
concludes that the RLF in high-z clusters has higher am - 
plitude than that of low- 2 clusters ()Branchesi et al. 2006) . 
Finally, a recent study shows there is evidence for a ~ 5 
times higher RAF in galaxies at z ~ 1 (M. Gladders 2005, 
private communication). Assuming this increase in RAF 
is representative of cluster galaxies, we adopt 7 = 2.5 for 
the power-law evolution model for the AGN. 

As an example of the degree of contamination due to 
AGNs, we show in Fig. ITU the cumulative fraction of clus- 
ters at z = 0.6 whose AGN flux is greater or equal to 
s\Ssze \ as a function of s, at 150 GHz. The chosen cluster 
mass is 2 x 10 14 M Q , close to the expected detection limit 
of some future SZE surveys (e.g. ACT, SPT). We consider 
the 7 = 2.5 case here. Two curves are shown; the dashed 
line is the estimated contamination fraction when the SID 
from all the sources for which we have flux measurements 
at both 1.4 and 4.85 GHz are used to extrapolate the RLF, 
and the dot-dashed curve is obtained when the SID of the 
BCGs is used (hereafter the all-SID and BCG-SID cases, 
respectively). In iJ2.4l we find evidence that the BCGs, 
on average, may have a steeper spectral index compared 
to non-BCG sources (the mean indices of the two popu- 
lations are —0.74 ± 0.05 and —0.47 ± 0.15, respectively; 
the distribution about the mean is broad for both popu- 
lations, however). Therefore, the extrapolated 150 GHz 
RLF based on the SID of the BCGs has a lower amplitude 
compared to the one constructed with the SID from all the 
sources, and results in a lower contamination fraction. For 
example, the LCFs (i.e. the value of the curves at s = 1) 
are 12% and 1.8% for the two cases. The Figure also sug- 
gests that about 32% (9%) of clusters host AGNs whose 
fluxes are at least 20% of that of the SZE, for the all-SID 
(BCG-SID) case. 

We do not attempt to model the difference in the SIDs 
in our Monte Carlo simulations. The differences between 
the dashed and dot-dashed curves in Fig. [JJ] represent the 
uncertainty due to our lack of knowledge in the spectral 



behavior of radio sources at high frequencies, with the 
BCG-SID case being an extreme. We expect a more so- 
phisticated treatment will yield estimates in the region 
bracketed by the curves. Below we will proceed with the 
RLFs constructed with the SID from all the sources, as 
they should be more representative of all the radio sources 
that may contaminate the SZE signal. 
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Fig. 15. — The 150 GHz lost cluster fraction (LCF), defined as the 
fraction of all clusters whose AGN flux is at least equal to their SZE 
flux. The circles, squares and triangles show the LCF at z = 0.1, 
0.6 and 1.1, as a function of cluster mass. The forecast is based on 
the simple power-law extrapolation from the 1.4 GHz model that 
we have adopted (see details in text). We assume a power-law evo- 
lution of AGN population with redshift. See text for discussion on 
the shape of the LCF. 

We note in passing that contamination due to super- 
luminous AGNs is very small (e.g. at 150 GHz, < 1% of 
all clusters are lost due to logP > 27 sources); the main 
sources of contamination are the AGNs at log P = 24 — 26. 

In Fig. [15] we show the LCF at 150 GHz as a function of 
cluster mass, estimated at three redshifts (circles: z = 0.1, 
squares: z — 0.6, triangles: z = 1.1). We only show the 
case for 7 = 2.5 (the no-evolution LCF is just a factor of 
(1 + z) 1 smaller at each redshift). At each redshift, the 
LCF first increases with cluster mass, reaches a maximum 
at masses 1 — 10 x 10 14 M Q , then decreases with mass. The 
latter behavior can be simply understood because, in our 
modeling, the total number of AGNs in a cluster simply 
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scales with mass, while the SZE flux has a steeper mass 
dependence (oc M 1 - 68 ); therefore at the high mass end the 
SZE flux wins over AGN. The decrease of the LCF toward 
the lower mass end, on the other hand, is caused by the 
shape of the AGN RLF. For simplicity, consider a cluster 
which hosts an AGN whose flux S is equal to the cluster 
SZE flux. The number of such sources per cluster is simply 
4>{P)V 1 where P is the luminosity corresponding to S, and 
V is the cluster volume. For a cluster ten times less mas- 
sive, its SZE flux is about 50 times weaker, and therefore 
one AGN with flux of 5*/50 can contaminate the SZE sig- 
nal. However, in the low mass system the number of such 
sources is (f>(P/50)V/10. For our extrapolated 150 GHz 
RLF, because of its mild slope (approximate power-law 
slope —0.5 ~ —0.6 for logP > 23), there are not enough 
AGNs in low mass clusters to contaminate the SZE to the 
same level as in the higher mass systems. 

From the figure one also notices that for high mass clus- 
ters, the LCF decreases with redshift, while the opposite 
happens for low mass ones; the LCF of intermediate mass 
clusters shows little redshift dependence. This also can be 
understood from the shape of the RLF. Consider a cluster 
at redshift z whose SZE flux is Ssze- At that redshift 
the AGN luminosity that corresponds to Ssze is P z - At a 
larger redshift z' , because the SZE flux is roughly constant 
with respect to redshift, the contamination radio luminos- 
ity is larger than P z by a factor of g (which is simply the 
ratio of the luminosity distances and fc-corrections at the 
two redshifts); at the same time, the amplitude of the RLF 
is increased by a factor of h = ((1 + z')/(\ + z)) 7 . There- 
fore, the relative size between the slope of the RLF and the 
ratio r — h/g would indicate the redshift behavior of the 
LCF. For low mass clusters, where the SZE flux is small, 
the portion of the RLF of interest is toward the faint end, 
whose slope is flatter than r, which means the increase 
of the AGN population with respect to z compensates 
the ever-increasing luminosity required to contaminate the 
SZE toward higher-z. For high mass clusters, because of 
the steepness of the RLF, the 7 = 2.5 evolution is not fast 
enough to maintain the level of contamination. Finally, for 
intermediate mass clusters, the balance between the shape 
of the RLF and our adopted AGN evolution produces the 
roughly redshift-independent LCF. 

In Fig. [TH we show the LCF for a 2 x 1O 14 M cluster 
as a function of redshift, at three frequencies of interest 
to some SZE experiments (30, 90 and 150 GHz; e.g. SZA, 
AMiBA, SPT, ACT). Only the 7 = 2.5 case is shown. It 
is apparent that the radio point source contamination is 
significantly reduced for experiments operating at higher 
frequencies. This is due to the increase of the SZE signal 
toward higher frequency, and also the tendency for cluster 
radio sources to have negative spectral indices. We note, 
however, in our estimates no attempt is made to model 
the removal of cluster radio point sources (cither through 
multi-frequency bolomctric observations or through inter- 
ferometric filtering) . This requires careful consideration of 
instrumental effects and observational strategies, which is 
beyond the scope of this paper. 

As we note earlier, two approximations have been used 
in our analysis (the extrapolation of the RLF based on the 
spectral behavior of radio sources at 1.4 to 4.85 GHz, and 
the power-law redshift evolution of radio-loud AGN abun- 



dance in clusters). Dedicated study of the spectral energy 
distribution of cluster radio sources, based on a large sam- 
ple of clusters that span a wide range in redshift, would be 
ideal to provide the data needed to carry out better fore- 
casts for the next-generation mm- wave SZE surveys. Such 
a knowledge is crucial, in that to extract precise cosmo- 
logical information from a survey, one needs to be able to 
control and model the systematics to few per cent level (e.g. 
iMaiumdar k Mohrl[200i iLima k Hull2005f ). Useful infor- 
mation about the model employed in the forecast may re- 
sult from a cross-correlation analysis of a well-defined clus- 
ter sample and WMAP maps (e.g. lAfshordi et all 120051 : 
iLieu et al.l[200il . 
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Fig. 16.— Lost cluster fraction for a 2 X 10 14 M Q cluster at 30, 90, 
and 150 GHz (from top to bottom), assuming a power-law redshift 
evolution of AGNs. For the simple power-law extrapolation model 
from 1.4 GHz that we have adopted (see details in text), at 150 
GHz the AGN contamination is ~ 10%, and the rate increases as 
the frequency of a SZE experiment decreases. 

8. SUMMARY 

We have studied the nature of radio-loud AGNs from 
a sample of 573 X-ray selected clusters, using the public 
data from the NVSS and 2MASS surveys. Without red- 
shift information for individual radio sources, we perform 
a statistical study of several ensemble properties of cluster 
radio sources, paying particular attention to any cluster 
mass-related trends. In particular, we investigate their 
radial distribution (in terms of the surface density), and 
the 1.4 GHz and i^T-band luminosity functions. By cross- 
correlating the NVSS and 2MASS catalogs, we study the 
fraction of galaxies that is active in the radio, and estimate 
the duty cycle of radio-loud AGNs. We also discuss the im- 
plications of our findings for the next-generation Sunyaev- 
Zel'dovich effect cluster surveys, as well as for the role of 
AGNs in the heating of the intracluster medium. Finally, 
we present evidence supporting a population of low mass 
(Mk > —23) cluster galaxies which host radio-loud AGNs 
(logP > 23, P in unit of WHz -1 ). In this paper we set 
an (arbitrarily chosen) division for high and low mass clus- 
ters at log Af2oo = 14.2. Our results can be summarized 
as follows: 
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1. We identify the brightest cluster galaxies (BCGs) in a 
statistical sense for a subsample of 342 clusters, based on 
the expected if-band magnitude and the location within 
clusters ( 323]) , We find that the 1.4 - 4.85 GHz spec- 
tral index distribution (SID) of the BCGs is different from 
that of non-BCG sources (i j2.4p ; in general the BCGs 
have a steeper spectrum. The mean spectral index for 
BCGs is a = -0.74 ± 0.05, while that for non-BCGs is 
a = -0.47±0.15. About 14% (17%) of BCGs (non-BCGs) 
show a rising spectrum (a > 0). 

2. The radial distribution of radio-loud AGNs is very cen- 
trally concentrated, and can be described by an NFW 
profile with concentration c ~ 25, when BCGs are ex- 
cluded. Inclusion of the BCGs increases the concentration 
to c > 50 ($3J Table [T|). More powerful radio sources 
(log P > 24) show a more concentrated distribution com- 
pared to the weaker sources ( N3.2|) . which may be caused 
by the luminosity seg regation of galaxies as s uggested by 
previous studies (e.g. iLedlow fe Oweniri995bD . or by the 
confining pressure due to the ICM. For Mr* (~ —24) 
or fainter galaxies, the distribution can be described by 
an NFW profile with c = 4 — 5; for the more luminous 
galaxies, however, the radial profile has c > 9 ( H3.3p . 
The concentration of galaxy distribution presented here 
is higher than that found in our previous study (LMS04); 
the main causes of the discrepancy are the cluster sample 
selection (only clusters not strongly affected by the back- 
ground large scale structure are used) and the cluster mass 
estimator. Our present analysis accounts for the variation 
in the background number count due to the large scale 
structure better. Finally, we find no indication for clus- 
ter mass-dependence of the radial distribution of the radio 
sources ( t)3.2p . 

3. After accounting for the difference in the mean over- 
densities in the field and in clusters (200/J7m), the am- 
plitude of the cluster radio luminosity function (RLF) is 
still higher than that of the (scaled) field RLF (gj}. Be- 
tween the luminosity interval logP = 24 — 27, the clus- 
ter AGN number densi ty is ~ 5700 times higher than 
the field value (see also iMassardi fc De Zottil I2004T ) , cor- 
responding to a 5700/(200/12m) ~ 6.8 times higher prob- 
ability for a galaxy to be radio-loud. An earlier study 
finds no difference between the field and cluster RLFs in 
terms of the fraction o f galaxies that is active in the radio 
(jLedlow fe Owenl fl996) ; in M4.4I we discuss possible origins 
of the different conclusions, and find that the comparison 
is non-trivial. 

4. The RLFs for high and low mass clusters show marked 
difference at the very bright end; the most luminous AGNs 
are found only in high mass clusters f ^4.1[) . The amplitude 
of the RLF for low mass clusters is slightly higher than that 
of high mass clusters; limited by the poor statistics, we do 
not further quantify the difference. 

5. We construct the if-band LFs for galaxies (based on 
2MASS sources) as well as for radio-loud AGNs (from the 
NVSS sources that have a 2MASS counterpart, which we 
refer to as the XC sources). By comparing the number 
of galaxies more luminous than a K-band magnitude for 
the 2MASS and XC sources (selected with a certain ra- 
dio luminosity range), we infer the radio active fraction 
(RAF) for cluster galaxies (JL2|). When the limiting radio 
(i-C-band) luminosity is increased (decreased), the RAF de- 



creases (see Table 2]). About 5% of cluster galaxies more 
massive than Mk* host AGNs. Our data also suggest that 
the RAF in massive clusters is a factor < 3 larger than that 
in low mass ones, with differences mainly coming from lu- 
minous galaxies (Mk < Mk*)- Because of the different 
radial distribution of radio sources and galaxies, we cau- 
tion that the RAF depends on the radius over which it is 
measured f £14.4|) . 

6. Special attention is paid to the RAF of BCGs ( jj^) . 
Compared to the majority of galaxies (e.g. Mk < Mk*), 
the BCGs are 10 times more likely to host an AGN of 
typical radio luminosity (logP > 23), and the probability 
enhancement of having a logP > 25 AGN is even higher. 
From another point of view, compared to the cluster galax- 
ies of comparable K-b&nd luminosity (but not centrally lo- 
cated in clusters), the BCGs still have 2 — 3.6 times larger 
probability of harboring AGNs. This implies that the cen- 
tral regions of clusters favor AGN activities. The BCG 
RAF shows a strong cluster mass dependence: 36% of the 
BCGs in high mass clusters have logP > 23, while only 
13% of their counterparts in low mass clusters are as pow- 
erful in the radio. 

7. Using the fraction of galaxies that host radio loud A GNs 
as a function of galaxy lu minosity (iBest et al.l |2005[) and 
the stellar mass function (|Bell et alj 2003f ) of field galaxy 
population, we compare the RAFs in clusters and in the 
field, and find that the latter is 10 — 20% of the former 

(E23). 

8. The RAF can be turned into an estimate of the typi- 
cal lifetime t r of radio sources ( q4.3[) . For Mk* galaxies 
formed at z = 3, t r w 5.7 x 10 8 yr, which is consistent with 
the quasar lifetime estimated from their spatial clustering 
(|Croom et alJ[2005h . 

9. By comparing the RLFs constructed from the XC 
sources and from the sources that are selected solely from 
the NVSS catalog, we infer that roughly 40 - 50% of the 
powerful cluster radio sources (logP > 23) reside in low 
mass galaxies (e.g. Mk > —23; see . Support for the ex- 
istence of this population comes from (i) their abundance 
as implied by the near-IR cluster LF, (ii) the surface den- 
sity profile and the RLF of the radio sources that do not 
have a 2MASS counterpart; in particular we show that 
they are strongly clustered. And finally (Hi) their optical 
color being consistent with the cluster red sequence. 

10. It is very likely that the central region of clusters plays 
an important role in fueling the central engine of radio loud 
AGNs; possible mechanisms include the higher confin- 
ing pressure of ICM, and the cooling instability ("cooling 
flow") towards cluster center. Observations that strongly 
support this view include (1) the concentrated spatial dis- 
tribution of radio sources, irrespective of the stellar mass 
of their host galaxies, (2) the RAF of BCGs is > 2 times 
higher than that of galaxies of comparable stellar mass, 
and (3) the RAF of cluster galaxies is larger than their 
counterparts in the field, and the amplitude of the clus- 
ter RLF is higher than that of the field RLF when the 
difference in the overdensities are accounted for. 

11. Because of their centrally concentrated spatial distri- 
bution, the radio-loud AGNs are an important heating 
source for the intracluster medium, especially within 10% 
of the virial radius. We estimate that AGNs can provide 
~ 0.13?7 keV per particle near cluster center, where 77 is a 



Radio Galaxies in Clusters 



23 



coefficient relating mechanical energy input to radio emis- 
sion luminosity (1/20 < i] < 20) (jg§. 
12. Using the SID derived at low frequencies (1.4 and 4.85 
GHz), we extrapolate our 1.4 GHz RLF to higher frequen- 
cies (^7]). Our extrapolation at 30 GHz agrees reasonably 
well with existing data. We provide fits to the RLFs at 
several frequencies in Table [5] Assuming a power-law evo- 
lution of the radio-loud AGN population, such that clus- 
ters of a given mass have five times more AGN at z = 1 
than at z = (but no change in the shape of the RLF), we 
model the contamination of cluster SZE signal by AGNs 
using a Monte Carlo scheme. We define the "lost cluster 
fraction" (LCF) as the fraction of clusters whose AGN flux 
is greater or equal to their SZE flux, Sagn > I-SszbI- At 
150 GHz, at a given redshift, the LCF first increases with 
mass, reaches a maximum of 10 — 13% at 1 — 10 x 10 14 M Q , 
then decreases with mass. For clusters more massive than 
10 15 M Q , the LCF decrease with redshift; low mass clus- 
ters (M 2 oo < 1O 14 M ) show the opposite behavior. For 
clusters of intermediate masses, the LCF decreases weakly 
with increasing redshift. We caution that our estimate is 
based on a large extrapolation in frequency for the spec- 
tral behavior of AGNs, and on the cluster AGN abundance 
evolution that is currently not well understood. Moreover, 
our estimates do not include the ability of cluster survey 
experiments to remove point source flux through angular 
or spectral techniques. 
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APPENDIX 

SYSTEMATICS AND UNCERTAINTIES IN SURFACE DENSITY PROFILE FITTING 

Radio Sources 

Here we provide some details about the fitting of the surface density profile of radio sources. First we argue for the use 
of a suitable cluster sample to study the stacked profile, then we describe the mock observations used to understand the 
systematics and uncertainties in the fitting procedure. 

In general, to study the density profile of sources statistically (i.e. without the aid of spectroscopic redshifts), one has 
to determine the background either globally (i.e. from the logA-logS") or locally (e.g. from an annulus surrounding the 
cluster where the contribution of cluster sources is negligible). At large radii, the estimates from both methods should 
agree, at least when many different line-of-sights are considered (that is, when many cluster fields are stacked). For 
example, for sources with log P > 23 in our RASS cluster sample, within r out = 2r2oo the locally derived background level 
is ~ 5% lower than that based on the log N-logS. With r ou t = 5r2oo, the two estimates agree to 0.2%. This suggests that 
the RASS sample is suitable for the study of the surface density profile of radio sources (this is not the case for 2MASS 
galaxies; see M3.3|) . 

We employ mock catalogs to help determine the uncertainty in the best fit profile (in particular that of the concentration, 
a c ) and understand the systematics of the fitting procedure. Given a cluster sample (e.g. the 188 clusters which allow the 
construction of surface density profile for sources with log P > 23 in Q and a concentration, for every cluster we generate 
cluster sources whose spatial distribution follows an NFW profile, in addition to uniformly distributed background sources. 
To incorporate the presence of large scale structure in the vicinity of cluster fields, for every real cluster we measure the 
surface densities of radio sources in annuli of radii 2 — 3r2oo and 3 — 5r2oo and compare them with the predictions based 
on the global number count, logiV-logS". Averaging over all clusters, the ratio of local number counts to the ones based 
on log AT-logS is almost unity in both annuli, with dispersion about 0.2, indicating that (1) the presence of the clusters 
do not affect the source count outside the virial region (an indication that radio sources are quite centrally distributed), 
(2) the background is fairly uniform and, on average, log N-logS gives good estimates of the background level (which 
confirms the conclusion drawn in the previous paragraph) . While the number of background sources in mock clusters is 
determined by the limiting radio luminosity, the log A^-log S 1 , and the angular extent of the cluster, we model the number 
of cluster radio sour ces N r by assum ing it is proportional to the number of galaxies N g in a cluster; using the -/V 9 -M 2 oo 
relation obtained byEffieLall (j2004j, LMS04), 

N g = 37(A/ 2 oo/10 14 M Q ) - 85 (for M K < -21 galaxies), (Al) 
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we write N r = 37^(M 2 oo/10 14 M Q ) - 85 . The actual number of sources assigned to a cluster is drawn from a Poisson 
random number whose mean is N r . We adjust the value of £ until the total numbers of sources within r 20 o and 5r 2 oo of 
the composite (stacked) cluster match well those in the observed composite cluster. The value of £ depends on the limiting 
radio power of sources; for logP > 23 sources we can construct realistic mock composite clusters by setting £ = 0.02. 
After choosing a suitable £, we generate 50 mock composite clusters and fit the resulting profiles; the standard deviation 
of the distribution of best-fit c is then used as <j c of the observed profile. 

We further use these mock observations to study the effects of uncertainties in the cluster mass (i.e. r 2 oo), cluster center, 
and position of individual radio sources on the surface density profile. For each cluster in the sample, we perturb its mass 
in logarithmic space from the one inferred from the X-ray luminosity by a Gaussian random number with a standard 
deviation of 0.2 (corresponding to ~ 50% fractional uncertainty in mass; q2.1j) . Because apertures of 3' r adius are used to 
locate the c entroid of extended sources from the raw RASS data in the NORAS and REFLEX surveys ()Bohringer et al.l 
2000. [2004T ). we expect the uncertainty in the center position to be < 1'. As our sources are stronger than 10 mJy, the 
uncertainty in the radio source position is about 1" fi j2.2p . We therefore randomly displace the cluster center and radio 
source positions up to 1' and 1", respectively. 

For objects with a very concentrated spatial distribution like radio sources, it is the uncertainty in the cluster center 
that has the largest impact on the shape (concentration) of the density profile. For a circularly symmetrical profile, false 
centers tend to smear the profile and lower the fitted concentration. The uncertainty in cluster mass translates into a 
scaling of r 2 oo; the concentration will be biased low (high) for clusters whose masses are underestimated (overestimated). 
Finally, as long as there is no systematic shift of the positions of the sources, the small positional uncertainty for NVSS 
sources should not bias the value of concentration. 

The effects of these uncertainties are quantified in Table [ATI We use 50 mock composite clusters, each is composed of 
188 clusters, to study how well our fitting code can recover the input concentration Q n in the Monte Carlo simulations, 
where we (1) assume there is no uncertainty at all (the "none" entry), (2) vary the cluster center up to 1' ("center"), (3) 
perturb the position of individual sources up to 1" ("individual"), (4) change the cluster mass in logarithmic space by 
a Gaussian with standard deviation of 0.2 ("mass"), and (5) combine all three uncertainties ("combined"). Three input 
values of concentration are considered (ci n = 20,25,30). We record the mean and standard deviation of the 50 best-fit 
values of concentration as cm in each of the entries. The number of mock observations is sufficiently large that the mean 
and scatter of cgt do not change much when we create more mock composite clusters. 

Table A7 

Effects of Uncertainties in Cluster Properties on Concentration 

(Radio Sources) 



Cfit 





none 


center 


individual 


mass 


combined 


30 


31 ±9 


28 ± 11 


31 ± 10 


32 ±21 


28 ±8 


25 


25 ± 10 


24 ±7 


24 ±6 


25 ±8 


24 ±6 


20 


21 ±6 


18 ±4 


20 ±5 


20 ±6 


19 ±5 



Note.- Results based on 50 mock composite clusters for sources more powerful than 
logP = 23. Ci n is the input concentration in the Monte Carlo simulations, c nt records 
the average and standard deviation based on 50 mock observations. See test for the 
meaning of different entries in the Table. 



From the Table we see that our fitting code is able to recover the input concentration without any biases, when no 
uncertainty is included in the simulations (the second column in the Table) . Very interestingly, the uncertainties in cluster 
mass only introduce a small, positive bias in the concentration. This is likely because cluster masses are perturbed high 
and low in the log space with similar probability; as long as the X-ray luminosity is a unbiased mass indicator over 
the mass range of clusters in our sample, such modeling of uncertainty in cluster mass is justified. More importantly, 
the effect of uncertainties in cluster center is slightly offset by the combined effect of uncertainties in cluster mass and 
position of radio sources, making the resulting profile only a little bit biased towards lower (4 — 7%) concentration. We 
note, however, if we perturb the cluster center up to 2', the bias on c will increase. Including only the central position 
uncertainty, on/cm ~ 0.7 for the three c m values considered. When all three kinds of uncertainties are present, the ratio 
is still about 0.7. Even though we could not rule out the possibility that some of the RASS clusters may have their center 
offset from the true location by as large as 2', the proportion of such clusters should be very small (< 2%; iBohringer et alj 
2000), and therefore our simulation would represent an extreme and very unlikely case. 

In summary, these tests suggest that the concentration of the spatial distribution of radio sources in clusters found in 
121 should be representative of the true value, with a possibility of being biased low by < 7%. 
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2MASS Galaxies 

We also use mock observations to study the effects of uncertainties in cluster properties (cluster center and mass) on 
the best-fit parameters of the galaxy surface density profile. The procedure is similar to that outlined in tjA.ll with some 
modifications that are described below. 

First, the number of cluster galaxies is based on Eq. IA11 N g = 37£'(M2oo/10 14 Mo) - 85 , where £' serves to scale the 
galaxy number to that appropriate for the chosen iGband magnitude limit Mkmiu- The i^-band luminosity function of 
cluster galaxies determined in LMS04 is used to set Second, the local-to-global number count ratios in the two annuli 
around each cluster (2 — 3r2oo and 3 — 5r2oo) show larger variation compared to the case of radio sources. Averaging 
over all clusters, the ratios in the inner and outer annuli are both greater than unity, with large dispersions. To better 
understand this, we have selected a subsample from all the RASS clusters which are not closer to each other than 6r2oo, 
and repeat the local-to-global number count comparison. With this subsample, the effects on the local background counts 
due to nearby clusters in projection should be much reduced (although there would still be clusters whose X-ray flux is 
below the limits of NORAS and REFLEX surveys present in the annulus regions). The ratios are 1.4±0.7 and 1.2±0.4 for 
the innter and outer annuli, respectively. In addition to variations in the large scale galaxy distribution, this also implies 
that the galaxies have a broad spatial distribution whose effect may be seen beyond 3r2oo! In fact, using a smaller sample 
of clusters which are isolated from each other by at least 8r2oo , we find that the local-to-global number count ratio in the 
annulus of radii 5 — 8r2oo is 1.0 ± 0.3. This exercise suggests that we should examine the spatial distribution of galaxies 
using data out to at least e.g. 5r2oo- To include field-to-field variations in the background galaxy counts, we incorporate 
the local-to-global number count ratio from the real clusters in the mock cluster observations. 

Combining our model of number of cluster galaxies, the logN-\ogS from the 2MASS all-sky data release, and the 
local background variation based on real clusters, we are able to generate mock composite clusters whose numbers of 
galaxies within r2oo and 5r2oo match those in the observed ones well. We further include uncertainties in the cluster 
center and mass in the mock clusters, by randomly displacing the cluster center up to 1', and perturbing the cluster 
mass in logarithmic space by a Gaussian random number with a standard deviation of 0.2. After creating an ensemble of 
composite clusters, we fit the surface density profiles and use the standard deviation of the resulting best-fit concentrations 
as the uncertainty of the concentration for the observed composite cluster. 



Table A8 

Effects of Uncertainties in Cluster Properties on Concentration 

(Galaxies) 









Cfit 






none 


center 


mass 


combined 


5.0 


5.0 ±0.5 


4.5 ±0.4 


5.0 ±0.6 


4.9 ±0.5 


4.0 


4.0 ±0.4 


3.6 ±0.3 


4.0 ±0.4 


3.9 ±0.4 


3.0 


3.0 ±0.3 


2.8 ±0.2 


3.0 ±0.3 


2.9 ±0.3 



Note.- Results based on 100 mock composite clusters for galaxies more luminous 
than Mk = —24.5. ci n is the input concentration in the Monte Carlo simulations, cm 
records the average and standard deviation based on 100 mock observations. See text 
for the meaning of different entries in the Table. 



We record in Table IA8I the effects of the uncertainties in cluster center and mass on the concentration. We construct 
100 mock composite clusters with galaxies more luminous than Mk = —24.5 (slightly more luminous than Mk*), with 
three values of input concentration (q n = 3, 4, 5). The second column shows the recovered concentration cgt, which is the 
mean and standard deviation of the distribution of the best-fit c from 100 mock observations, when no uncertainties are 
included. The third and fourth columns are the results when the cluster center and mass are perturbed, respectively, and 
the last column shows the result when both uncertainties are present. Our code can recover the true concentration without 
any bias when the mass and cluster center are perfectly known. Any error in cluster center determination translates into 
a bias in the concentration. The uncertainty in cluster mass, on the other hand, does not produce any bias, as discussed 
in £|A.ll This test shows that the concentration we measure would be very close to the true value, with a potential bias 
of < 4% toward lower value. 
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